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Interesting websites: 

The following are two excellent websites from German teachers using Derive in Math 
education. On Mr. Grude's homepage you can find a very fine contribution on Cellular 
Automata. 

http://home.arcor.de/gallinat.richard/arcor/index.htm 

http://www.dietrichgrude.de/kranich/indexkranich.htm 


Our French friends produced an online course for Tl-calculators: 

http://www.t3ww.org/france/t3-online.htm 

Gerhard Hainscho produced a very extended Voyage200 Course for Beginners. Walter 
Wegscheider put this course as an online course on the web (together with many additional 
resources and links). Many thanks Gerhard and Walter from T -Austria and ACDCA. 

http://www.austromath.at/daten/voyage200 
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p i 


Dear DUG Members, 


this DNL does not contain so many but more extended contributions. Revising DNL#3 and preparing 
its pdf-publication I got the idea to treat the Differential Equations from 1991 with the TI- 
92+/Voyage200. I gave a workshop on DEs based on course materials provided by Gunter Redl, so I 
added this paper too (together with some English explanations). I hope that the DERIVE users are not 
too disappointed, but we didn't have so many CAS-Calculator contributions in the last DNLs. 

Many tasks in the WS are provided for students and WS-participants. The 2nd part (page 20) refers to 
DNLs#2/#3. You can download both in revised and updated form as pdf-files and there you will find 
more explanations. It was fascinating to apply the DE-Mode of the TIs , to compare the results and in 
case if the TI's DE-mode fails to program the respective functions according to the auxiliary functions 
provided in Derive 's ODE1.MTH. 

The second extended article is the result of a Swedish-Austrian cooperation between David Sjostrand 
an me. Last summer David presented his Excel/Derive cooperation on Cellular Automata in Gothen- 
burg. This was a challenge for me to wrap the procedures in a program. Finally we decided to give a 
joint presentation at the Matematik-Bienallen in Malmo, last January. We had much fun in the audi- 
ence and we would like to share our product with you. I'd like to call your attention to D Grude's 
website (see Information). 

On page 43 a tool is offered which might overcome compatibility problems between DERIVE6- and 
DER1VE5 files. I am very grateful for many valueable co mm ents and support from Theresa Shelby and 
Albert Rich. 


Finally after a long time we have again an ACDC contribution (Amazing Comer of the Deriver's Cu- 
riosity, founded by Alfonso Poblacion in DNL#22, 1996), which was sent by Milton Lesmes from far 
Colombia. Muchas Gracias, Milton. 

Many new members have joint the DUG during the first months of 2004. A warm welcome to you all. 
Torbjom Aim sent a contribution about Pollard-Factorization, which will be published soon. Among 
others I received an article from H.R. Geyer, who presents a roboter in 3D using the slider bars to 
move his "little man", and an extended article dealing with "Actuarial Mathematic" from MacDonald 
Phillips. 

In the meanwhile an interesting discussion was raised about the function value of sign(0), which is not 
finished. I'll report about it next time. 


I'll close repeating the invitation to join the DUG meeting in Montreal and wishing you all a wonder- 
ful summer. I am looking forward to meeting you with DNL#55 in September. 



DUG-Meeting 2004 
TIME-2004, Montreal 
Saturday, 17 July 2004, Lunchtime 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 77 User Group. 
It is published at least four times a year 
with a contents of 44 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE 
and the TI-89 1921 Voyage 200 as well as to 
create a group to discuss the possibilities 
of new methodical and didactical manners 
in teaching mathematics. 

As many of the DERIVE Users are also 
using the CAS -TIs the DNL tries to com- 
bine the applications of these modem tech- 
nologies. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 


Next issue: September 2004 

Deadline 15 August 2004 


Preview: Contributions waiting to be published 

Finite continued fractions St. Welke, GER 

Some simulations of Random Experiments, J. Bohm, AUT 

Wonderful World of Pedal Curves, J. Bohm 

Another Task for End Examination, J. Lechner, AUT 

Tools for 3D-Problems, P. Luke-Rosendahl, GER 

ANOVA with DERIVE & Tl, M. R. Phillips, USA 

Hill-Encription, J. Bohm 

CAD-Design with DERIVE and the Tl, J. Bohm 

Avoiding Convolution and Transforming Methods, M. Lesmes-Acosta, COL 
Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Fuzzy Logic, G. Hagen, AUT 

Simulating a Graphing Calculator in DERIVE, J. Bohm, AUT 
Modelling real date: Enthalpie Values 
A TVM-Calculator for DERIVE, M. R. Phillips, USA 
Shares, Put and Call, J. Bohm, AUT 
Actuarial Mathematics, M. R. Phillips, USA 
A Roboterwith Derive, H.R. Geyer, GER 

and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 

URLs: http : //www . austromath .at/ dug/ index, htm 

http : //shop . bk- teachware . com/k . asp?session=2579115&kat=ll 
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Adam Marlewski algebraistxxi0yahoo.com 

Hello All, 


Lou Lowell wrote: 

I recently examined the graph of the expression (x A 2+5x+6)/(x+2) in both Dfw4.11 and Dfw 5.06. In 
ver 4.1 1 the discontinuity at x= -2 is observable on the graph, however the same function in ver.5 is 
never observable. Moreover, in ver.4 when tracing along the curve, derive gives a little click when 
passing the discontinuity; in ver.5 there is no indication of a discontinuity at -2. 

Can this be changed to be visible in ver.5? 


I answer: In DfW 5.04 the discontinuity at x = -2 is visible if the scale is not too big. See atached file 
(_discont.dfw) please. 

Adam 


DNL: I'd like to add another comment: It makes a difference if you have set Option Simplify/App- 
roximate before plotting in the 2D-Plot Window or not. If yes, Derive automatically simplifies the 
x 1 + 5x + 6 


expression 


x + 2 


to x + 3 and subsequently plots the line without any gap. 




MacDonald Phillips phillipsm0GAO.GOV 

The problem submitted by Jan Vermeylen on page 18 of DNL#53 presents quite a problem in equati- 
on fitting. While the resulting, very bulky, equation is linear in k, it is definitely not linear in K. Since 
there are only two parameters to estimate, it is possible to use any two rows of data and set up two 
equations in two unknows to obtain initial guesses for the parmeters. However, in solving the equati- 
ons for the initial guesses, k and K turn out to be complex values. This presented a problem since my 
nonliear regression routines in the User directory of Derive 6 did not handle datasets with complex 
numbers or complex parameters. So I modified them so they do! I'm sure complex answers were not 
expected, but given that there are only three data points, that's the best I can come up with at this 
time. 

So, I get the following: k = 0. 041708+0. 5551 27*i, K = -43. 2998+45. 649*i where i is the square root 
of -1 . 

The adjusted R A 2 is only .75 and the t and F statistics indicate that the fit is not statistically significant. 
But that's to be expected with only 3 data points. 

Any other thoughts on solving this problem would be greatly appreciated. 

If anyone wants a copy of my latest regression routines (they also handle weighted regressions now), 
please email me at phillipsm@gao.gov. 


Don Phillips 
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Albert Rich 

Hello Josef, 


Theresa ask me to respond to the question raised in your email that follows. In Derive 5, when as- 
signments are simplified, the result is just the simplified right side of the assignment. In Derive 6, 
when assignments are simplified, the result is the assignment with the right side simplified (this is 
stated in the Derive 6 on-line help for the Author > Variable Value command). For example, when the 
assignment 

x := 2+3 


is simplified in Derive 5, 5 is returned; whereas in Derive 6, the assignment x := 5 is returned. This 
change in the behavior of Derive was made for two reasons: 

1. When an expression is simplified in Derive, the result should be equivalent to the origi- 
nal expression. In the example above, 5 is not equivalent to the assignment x := 2+3; whereas, x := 5 
is equivalent to it. 

2. In Derive 5 if an assignment is entered and simplified on the Expression Entry line, only the simpli- 
fied right side of the assignment appears as an expression in the algebra window. In Derive 6, the 
entire assignment appears as an expression in the algebra window. Thus if the algebra window is 
saved as a dfw or mth file, it is clear that an assignment was made when the file is subse- 
quently loaded. 

However, as you discovered, for efficiency reasons, when assignments to local variables made within 
function definitions are simplified, only the simplified right side is returned. 


Hope this explanation helps and I apologize for the inconvenience. In the next release of Derive 6, 
I will clearly state this change in behavior of Derive 6 in the on-line help. 

Aloha, 

Albert 


Tom Miles temiles0cox.net 

Sharing some interesting observations 

Hello Josef, 

I thought I would pass along some interesting observations I have made regarding the solutions to my 
equations. 

As I mentioned in a previous communication, I was frustrated that my programs were failing when I 
was trying to solve for polygons with sides greater than 40. I had been using Derive 6 in the decimal 
mode. To reach a solution, I often had to increase the decimal precision to over 100. I made two 
changes that have greatly improved the situation: 

First, I converted all the equations to the radians and changed the system to Angle := Radian. 

Second, I simply changed NSOLVE((y - yl) = m(x - xl)) to NSOLVE((y - yl) - m(x - xl) = 0). Sur- 
prisingly, setting the relationship to “0” made a big difference. 

Hope this information may be of interest. 

Sincerely, 

Tom Miles 







Valeriu Anisiu 

Hello Albert, 

I have some Derive problems for you. Some of them were not present in version 6.01 Beta. 


l 


7(2 

0 

For the antideriavative a strange branch is chosen. Same problem with the next integral: 
l 

1 73-LN(4-76 - 5-./3 - 6-72 + 10) n 

2 dx = - - — + 72 

7(1 - x) + 7(x + 1) + 1 6 4 

0 

3 After a long time a complicated multiple SUBST ex- 
pression is returned; of course the integral can be ea 
sily computed (e.g with INTSUBST). Derive should 
either quit (as in other cases) or give a normal result. 

4 Memory exhaustion in either exact/approx mode. Also for other sums related to ZETA(). 

r 2 

73-LNCt-(72 - 1) - 73 + 72) 73-LN(t-(72 - 1) + 73 + 7 2) 7(2 - x ) - 72 

c- SUBST ATAN(t) + , t, 

° L 6 6 x 

exhaust memory! 

6. (Interface problem) 

In Derive 6.0 (and 6.1 beta) the mouse wheel did not work; now it works but in the upward direction 
the worksheet jumps very fast to the top (being almost useless). 

I finish with a question related to the approx mode in Derive. 

From your very clear explanation I understood that Derive does not use a floating point arithmetic. 

A real number seems to be approximated with PrecisionDigits=n as A/B*10 A C with abs(A)<K1(n), 
0<B<K2(n) and no limit for the integer C. [possibly 2 A C is used instead of 1 0 A C]. 

I am curious to know: 

- which are the values for K1(n) and K2(n)? Are not they too "optimistic"? 

- why a limit (maybe optional) for C is absent (for overflow/underflow in approx mode)?. 

Because of this, in numerical computations "infinitesimals" may propagate. Derive behaves some- 
times strangely when abs(C) is big. For example ITERATE(f(x),...) may appear to diverge even if the 
iterates of f converge to 0 or to a vector having a 0 component. 

I wish you a happy Easter [this year the Western and Eastern Easter coincide as a Derive program 
posted a couple of years ago shows]. 


1 

dx 

2 

x • 7( x +x+l)+2-x+2 




x 73-LN(4-./6 - 5-73 - 6-72 + 10) n .[2 

dx = - - — + - 

2 12 826 

- x ) + x + 1 


Valeriu 




D-N-L#54 


Derive- and CAS-TI-User Forum 


Albert Rich 

Hello Valeriu, 

The following is in response to your numbered questions: 

1. Unfortunately, as plotting both the real and imaginary parts shows, the antiderivative Derive finds 
for x/(SQRT(2-x A 2)+x+1) is discontinuous at x=SQRT(3)/2-1/2. The correct answer to the definite 
integral is 

SQRT(3)*LN(-4*SQRT(6)+5*SQRT(3)-6*SQRT(2)+10)/12-pi/8+SQRT(2)/2 

This result can be obtained by adding the integral from x=0 to the discontinuity to the integral from the 
discontinuity to 1 . I will add this problem to my list of discontinuous antiderivatives to resolve. 

2. One application of Display Step shows that the antiderivative of 1/(SQRT(1-x)+SQRT(x+1)+1) re- 
duces to the problem described in #1 above. 

3. Unfortunately, as the Display Step feature shows, Derive rationalizes the denominator and then 
applies partial fraction expansion resulting in explosive growth of the problem. Can you tell me what 
substitution is required to find the antiderivative? 

4. In the next release of Derive this summation problem will almost instantly simplify to pi A 2/6 - 
ZETA(2, 50001) and approximate to 1.644914067. 

5. This memory exhaustion is due to a stack overflow problem that occurs when trivial factoring some 
algebraic expressions. For example, try trivially factoring x*(2*SQRT(6)+2*SQRT(2)-8)+4. I hope to 
have this resolved for the next release of Derive. 

6. On my computer running Derive 6.01, scrolling the algebra window using the mouse wheel works 
normally. I will ask Theresa if she can reproduce and resolve the problem you are having. 

Concerning your questions concerning rational approximation in Derive: 

Derive uses the approximate rational arithemetic provided by muLISP (the LISP language pseudo- 
code interpreter in which Derive is written). muLISP approximates a number by finding the clos- 
est rational approximation for which the sum of the bit lengths of the numerator and denominator does 
not exceed a given number of bits. For n digits of precision set by the user, Derive uses an internal 
precision of 4*CEILING(5*n+6,6) bits. Thus for the default 10 digits of precision, the sum of the bit 
lengths of the numerator and denominator is limited to about 40 bits. However, there is no limit (other 
than memory) on the bit lengths allowed for integers and pure reciprocals, since this is the only way to 
preserve their magnitude. 

As you noted, there is no underflow or overflow in Derive other than memory exhaustion. Overflowing 
to infinity is a bad idea for a computer algebra system, since even big numbers are much closer to 1 
than they are infinity. A golden rule (or at least ideal) in Derive is that 1/(1/n)) should equal n. If for 
n»0, 1/n approximated to 0, 1/(1/n) would approximate to infinity. However, as I say above even if 
n»0, n should not approximate to infinity. Therefore, it is also a bad idea to underflow to zero. 

Hope this explanation helps. I have copied this email to several others since the discussion of ap- 
proximate rational arithmetic may be of interest to them. I hope you do not object. 


Aloha, 

Albert 
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Valeriu Anisiu 


Hello Albert, 

here are some remarks related to your last email. 


INT 


For an integral of a function R(x,sqrt(ax A 2+bx+c)) where R is a rational function, it is not a good idea 
to rationalize the denominator; the Euler's substitutions are recommended. They are implemented (as 
a bonus) in INTSUBST (because I noticed that they are not implemented in Derive), 
e.g. 

INTSUBST(1/(x*SQRT(x A 2+x+1 )+2*x+2),x, x A 2+x+1="euler1") 

— > 

-SQRT(3*SQRT(1 7)/1 7+5/1 7)*ATAN((4*SQRT(x A 2+x+ 1 )-4*x-SQRT(1 7)+3)*SQRT(1 1 *SQRT(~ 
17)/64+45/64))+SQRT(3*SQRT(17)/68-5/68)*LN(4*SQRT(x A 2+x+1 )-4*x-SQRT(22*SQRT(1~ 
7)+90)+SQRT(17)+3)-SQRT(3*SQRT(17)/68-5/68)*LN(4*SQRT(x A 2+x+1 )-4*x+SQRT(22*SQ~ 

RT(1 7)+90)+SQRT(1 7)+3) 

(the change of variable is sqrt(x A 2+x+1)=x+t) 

APPROX 


I don't think that for numerical computations, over/underflowing is such a bad idea (of course I do not 
suggest to change Derive for this but it is anyway largely used). Your statement that a huge number is 
closer to 1 than to inf is also arguable because in the interval [-inf, inf] usually other metrics are used 
(for example via a bijection to [-1,1]). 

The fact that in approx mode 1/(1/n) = n is of course fine but why should it be true only when n is inte- 
ger; in approx mode this is not essential for a user who wants seriously to approximate and not to 
"test" 

superficially Derive. 

I have recently seen a proposal for approximating real numbers in rational arithmetic where x is ap- 
proximated to the first convergent p/q (of the continuous fraction of x) for which abs(x- 
p/q)<=eps*abs(x). 

With this approach PrecisionDigits would have a full meaning; the problem of the uneven distribution 
of Q is resolved with a very moderate cost! 

The kernel functions should return the result within the requested precision. In Derive (unlike any 
other CAS I know) this is not the case with trigonometric functions. E.g.: 

APPROX(sin(10 A 11),10) — > 0.7071067811 (which seems to be sqrt(2)/2 ??) instead of 
0.9286936604 

In this case the argument of the sine must be reduced internally using PrecisionDigits+1 1 
A very rough user solution could be to define: 

sin_(x):=IF(IDENTICAL?(Precision, Exact), 

SIN(x), 

SIN(APPROX(MOD(x,2*pi),PrecisionDigits+MAX(LOG(ABS(x),10),0)))) 


Salut, 

Valeriu 


You can find the end of this discussion on page 45. 






D-N-L#54 


A Differential Equations Workshop for Handheld 


P 9 


Differentialgleichungen mit dem TI-89/TI-92+/Voyage 200 

Differential Equations with TI-89/TI-92+ and V200 

A Workshop 


Hier soil nicht die Theorie der DGL und ihre Losungsmechanismen besprochen, sondem an einigen 
ausgewahlten Aufgaben die Vorgangsweise auf den TI-CAS-Rechnem beschrieben werden. Der 
TI-89/TI-92PLUS und der Voyage 200 verfugen iiber ein gehoriges Potential an graphischen und ana- 
lytischen Werkzeugen zur Behandlung der Differentialgleichungen. 

We will not deal with theory of DEs, but show the solving procedure using the TI-CAS calculators. 


Beispiel 1/Example 1 

y=v(x+i) 

a) Zeichne das Richtungsfeld. Plot the direction field. 

b) Zeichne spezielle Losungen der DGL fur y(0) = 3, y(0) = -3 und y(-l) = 4 

c) Bestimme die allgemeine Losung und die oben angefuhrten speziellen Losungen mit Hilfe des 
Rechners indem du die handische Vorgangsweise nachvollziehst. 

Find the general solution reproducing the manual solving procedure. 

d) Bestimme diese Losung mit alien vom CAS zur Verfugung gestellten Hilfsmitteln. 
Die Losung aus b) konnte mit der entsprechenden nun gefundenen Losung iiberlagert werden. 
Use all means of CAS to find the solutions. 


a) Der Rechner wird im MODE auf Graph DIFF EQUATIONS eingestellt. Die Differentialgleichung wird in 
den [Y=]-Lditor geschrieben. Aber Achtung: die Rolle des x wird vom Parameter t ubernommen! (x 
must be replaced by f!) 



m 

zSSJ 1 

rfiT^DT y f3 y fh y fet 

| ▼ 4 — |Zoon|Trace|Regraph |Math 

raviic i 

J-LDTS 

iO=G. 

^yl '=yl (t + 1) 
yil=l 

y2' = 

yi2= 
y3' = 
yi3= 

y4' = 

yi4= 
y5' = 

LG=ai 
tnax=3. 
tstep=.G5 
tplot=0. 
xnin= _ 4. 
xnax=2. 
xscl=l . 
ynin= _ lG. 
ynax=10. 
yscl=l. 
ncurves=0 . 
Esteps 1 . 

y \ \ \ \ \ \ 

\ \\\\\\ v ^ ^ 

\ \ \ \ \ \ "v - 

\ \ \ Y W--------- 

/'///iff 
f / / / f / f 
'////// 
-- s' s' S S S / 

/ / X s' s' s' ^ --- ■ 

/ / / / / / S s ^ ^ 

f / / / / / / S ^ ■- *- ^ ^ 
iff//// w \ 

-- 'v. \ V \ \ 

\ \ \ \ \ ) 
WWW 1 

\ \ \ \ y \ 

yil= 

SELECT DUE 1ST-0F:DEFI FUNCTION ONLY 

MHIN 

F;HD HUTU DE 

MHIN Fifi[i HUTD DE 


Wenn Du in den [WINDOW]-Werten fur ncurves etwa den Wert 6 einsetzt, dann werden gleichmaBig 
iiber den y-Bcrcich 6 Losungskurven eingezeichnet. fiber F3 Trace konnen die Kurven abgetastet 
werden. 



b) Um die beiden Losungskurven zu den gegebenen Anfangsbedingungen zu sehen, gibt man die 
beiden Funktionswerte zu x = t = 0 in einer Liste zu yil an. Fur die dritte gesuchte Losungskurve 
muss der Wert fur tO auf-1 und yil auf 4 gesetzt werden. 


p 
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s \ \ W\1 
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. . . . \ i \ 

MAIN FiHti HUTD PE 

USE Y+U UR TV PE + [ESC]=CANCEL 

MAIN RAD HUTD DE 


Es gibt noch eine weitere - interaktive - Moglichkeit, Losungskurven durch bestimmte Punkte zu 
legen. Uber F8 1C konnen Anfangsbedingungen (Initial Conditions) direkt eingegeben werden. So 
wollen wir zB die Kurve durch P(-2|-2) legen. 


l?i2feS n T T ^c,fe.<?«*llShlDK' U F2 \ 

l?i2fcS’ n Ti^fe.<?«*hShlDK' u F2 SW 1 

y \ \ \ V \\ A 

u \ \ \ \V ^ ^ ^ ^ ^ ^ 

\ \ \ A A A 

A v N \ w--------- --- . 

V / / / t i 
r/////\ 
r / / / / / / 

Z Z Z / 

y \ \ \ \ \\ A ^ 

) \ \ \ \ \\% ^ ^ ^ / 

\ \ \ A A A 

A x x N — --- 

■ 7 / / / i 
r/ / / / / i 
r / / s / / / 

- z s s / 

/ / / / / / s -- ----- 

initial tonditioris=?' ^ ' v 
t=-2 . . yl=-3 

x -v A 
^ A A A \ \ J. 

\ A \ \ \ \ i 

. . . . \ i \ 

/ / / / — ^ 

l / / // f / / s ^ -* -- ^ 'v 
iff// / / / / A 

■%, 'x % A 

k \\\ \ \ \ 

k A \\ \ \ i 
; \ \ \\\ i \ 

MAIN RAD HUTD DE 

MAIN RAD HUTD DE 


c) Als Losungsmethode kommt hier die Methode der Trennung der Variablen in Frage. 
We apply Separation of Variables. 

— = y(x + 1) — > — = (x + 1) dx 
dx y 


upT 1 

u P n 

If- Ifl 1 qebra fcaTclotherlpr-qn I O Ic lean UpH 

■ f[T-)ay =tc* + l)ax ln(|y|)=-^- + x + c. 

f x 2 1 

■ sol vel ln(y) = + x + c , yl 

x 2 

-R- + X + c 

y = ■? *■ 

■y = f 2 ‘ 1 x = 0 and y = i -3 3J 

C -3 33 = <? c 

x 2 

9 + X 

■ ~Z e £ -^ylCx) Done 

X 2 

■3 ■€■ * -»-y2Cx) Done 

■c-Solue({-3 3J = 
o= Cln(3) + (2 fin: 

x 2 

J^ + x + c 
■ y = e z 1 c 

-1 

= f c , c) 

2-n + iOl 1 n(!3) + 2 ■ En3 
= 0 lnC3) + (2 ■ Gn2 ■ k + x) ■ * 

r x 2 ^ x 2 i 

l "3 ■ y” 2 ” 3 ,~ A 

solue <ln<y>=x >v 2/2+x+c,y> 


3*e A <x A 2/2+x>4v2<x> 

|y=e A (x A 2/2+x+c>lans(i> 


SIMLLATI RAD HUTD FUNC E/30 

1 SIMULATI RAD HUTD 

FUNC !Y30 


Alle drei speziellen Losungen werden gefunden und uber den Funktioneneditor gezeichnet. 
Beachte den Einsatz der cSolve-Funktion (wenn man nicht direkt fur e substituiert). 

We find all special solutions and plot them. Notice application of cSolve if not substituting directly 



d) Mit der implementierten deSolve-Funktion lasst sich die DGL direkt losen. Auch hier fuhrt cSolve 
zum gewiinschten Ergebnis. Die Zeichen stehen fur freie Parameter. 

The @- characters are free parameters. 


| - |r l qebrabaT:- [other Ipnqn I ole laan UpI 1 

f- |r i qebralcaT:- [other Ipnqn I oElean UpI 1 

I- Ifl 1 gebralcaTialotharlpnyn 1 0 Ic lean UpI 1 

■deSolve(y’ = y -(x + 1) , x , y) 

x 2 * 

~9~ + x 

y = Gl-tf- 2 

f TT + X I 1 

■ solue| y = i?l e *- | x = - 1 and y = 4 ; GlJ 

ei = 4 -li 

■ solvef 

■ cSolve 

X 2 

— + x ► 

y = @l-e ^ | x = 0 and y = t -3 

false 

f 4 + *, ► 

[y = G1 * | x = G and y = € -3 

G 1 = C "3 3J 

■deSolveCy 1 =y(x+l) and yCG)= - 3,x,y) 

x 2 + 
-=- + x 
y = -3 • * ^ 

■deSolvefy' =y-(x + l) and y(G) = 3,x,y) 

x 2 + 
~9~ + X 

y = 3- e- ^ 

solue (ans <1 > 1 x=“l and y=4,Pl> 

. A 2/2 + x> 1 x=0 and y=<-3 .3>,P1>| 

,,,y*Cx+l> and y<G>=3,x,y>| 

SIMULATI RAD AUTD FUNC 

SIMULATI 

RAD AUTD FUNC lEf'3D 

MAIN RAD HUTD FUNC il^O 


Auch die Randbedingungen (oder andere gegebene Punkte der Losungskurve) konnen sofort mit 
angegeben werden. We can include points of the integral curve into the deSolve-function. 
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Beispiel 2/Example 2 

/ = 

x V x ^ 

a) Bestimme die Losung mit Hilfe des Rechners, ohne vorerst desolve() einzusetzen. 

Tipp: Substituiere z = — . Solve without using deSolve() - substitution!!. 

x 

b) Wie lautet die spezielle Losung fur v(l) = %[e ? Special solution for y(1) = \fe 

c) Zeichne das Richtungsfeld und lege die gefundene Losungskurve driiber. Direction field. 

d) Bestimme diese Losung mit alien vom CAS zur Verftigung gestellten Hilfsmitteln. 

e) Fixiere die Losungen von a), b) und d) in einem Script (Textfile). Produce a script. 


a) z = - — > y - z-x — > v' = z'-x + z 
x 

Daher wird aus der DGL: z' ■ x + z - 2z ■ ln(3 z) —>—■ x = 2 z • ln(3z) - z — > — = — . 

dx 2 z ■ ln(3z) - z x 

Die Trennung der Variablen ist gegliickt. Es ffagt sich nur, ob die Integrate geschlossen ausgewer- 
tet werden konnen? Separation was possible, is integration possible, too? 




^ ^ * J > - l-F'l-.l'l ■ - 


iva . ? VVl A-al ,m ; I 






'.f[ 2 ■z-lrO z)-z h=J[|]^ 

1^2 W3 zi-0 =ln(|x|) + |: 
, 1 n( 2 • 1 n( 5 • z) - 1 ) 


= ln(|x|) + c 


- = ln(|x|) + c, y 


- 1 x = 1 and y = *■ ■ * 


- = ln(|x|) + c| z = — 


,3*g>~l>/2=ln<ahs<x>>+c I z=v/x 


FifiU BUTU FUNC H,'3P 


FifiU HUTD FUNC H . 1 B 


Das Integral lasst sich auswerten. Welche Integrationsregel(n) fuhren zum Ziel? [1] . 
Which integration rules are leading us to success? [page 121 


1 t I fl 1 qebralcalclother Iprffn I olcieall UpI 1 



— + lz2 

■solve ^ ,c 

In 
c = — 

^ 2 ‘ C -X 2 

f 61n(3)- n 

3 |c= 2 

^.(1^35 - l,'6)x 2 + 1x2 
y 3 

x .JInCS)-l/6V>c 2 + U2 

■ 1 ^ + yltx) Done 


J 

3 J 

2 

6- lnf3^ - 1 1 

xc: 1 . 

yc: 1.3956 



MHIN F;HD HUTD FUNC H.-'SO 

MHIN FifiU HUTD FUNC H/30 

MHIN F:HD HUTD FUNC 


b) 


c) 


Es gibt eine allgemeine und die gesuchte spezielle Losung. Der 
Graph der Losungskurve wird als Bild FI 2:Save Copy As unter 
dem Namen beisp2 abgespeichert - um ihn spater liber das 
Richtungsfeld zu legen. We save the graph for later combining 
it with the direction field. 









■■■ 

SAVE CDPV AS A 



Type: 

Folder: 

Uariable 

Picture-* 
main-* 
IbeispS II 





CEnt-er=SRUE > 

CESC=CRNCEL"> 



/ 



MHIN 

FiFlD HUTU 

FUNC 


fiber MODE wird in den DIFF EQUATIONS Modus gewechselt. Die [WINDOW]-Werte sollen mit denen 
von vorhin identisch sein. 


Pi=tss’ n TE?itr'Mnfeffi.fc^..i l 

I^Iz^T^cAe^aPhlM^hlD^XVllC I 

I^Iz^TricelRe^aPhlnJthWVllC | 

JLDTS 

to=i. 

y . , 2-yl . f 3 yl 1 

yl = t ln l. i J 

yil=l 
y2' = 
yi2= 

y?;= 

yy3= 

y4' = 

y 

III III!!// 
Ill III//// 

1 1 S It//// / 
1 ./ / / //- 

y 

III//////// 
III III//// 

1 1 S 11///// 

1 t fjf///'- — ---- 

/ / 1 1 

//////II II 
////III 1 1 1 
/////// 1 1 1 

t 

'/’■/j^Ti i i 

//////it ii 

/ / / t Si! /III 
//III!/ 1 1 1 

t 

yii= 

SELECT ONE iST-URDEF FUNCTIUN UNLV 

MHIN FifiU HUTU IE 

MHIN FifiU HUTU IE 


fiber FI LOpen wird beisp2 fiber das Richtungsfeld gelegt - und alles passt. 
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d) Wenn es "nur" um die Losung geht, versucht man gleich desolveQ. 


[^iHlSraKtUh^lp^IoJci^ UpI 1 

[^iHlSratftb&Ip^IohaST UpI 1 

(InO- 1 / 6 ) x 2 + 1/2 

■ ^ + yl(x) Done 

■ deSolueJy' = ■ 1 , x , yj 

y = 2 ^ and x > 0 

, 2 'H.T 

2 +UZ ► 
■ y = -^-1 3 and x i 0 | x = 1 ai 


C1 2 

solu6(ans(i),y> 

ans<l>lx=l and y=e A <i/3>| 

MAIN FsrtD HUTU DE 9230 

MAIN RAD HUTD DE 11/30 


U pf i 




3 l el_ 2 

„Cln£S - 126 ) x 2 + 122 


... e (2*Pl>*x 2/2+l/2>/3 I ans(i> 


DE 12^30 


e) Der komplette Prozess wird in einer Textdatei gespeichert (FI 2:Save Copy As). 





rFiT^Y f2t Y F2w T FH T FE T 1 

▼ f — | C o m m a n d | U i a w | E x a c u t a |F i nd . . . | | 

I'Fi-nriY F2T Y F3,F Y fh T fe Y i 

I t * — |uon man d |U i aw lExecute |F i nd... | I 



f OPEN ^ 


C:K<z{x)*x, x? 

C: S< lzc:2*z*ln<3*zI>-zD „ z)=m/x,x, cD 
C: ln<2*ln{3*z>-i::'---2=ln<abs<x>>+clz=yYx 
C : so 1 y e <1 n C 2* 1 n < 3*y zx > - 1 V2= 1 n < abs < x ) ) +c 
, y) 

C:y=x>t<#- A <:e A <:2^c.^x A 2z2+lz2V3l x=l and y= 
a A Clz3> 

C:solyeC^ A az35=^C^C2*cV2+lz2V3 J c> 

C : y=xYt-" ■:] (2*c) *x A 2/'2+ 1 / 2 V3 1 c= 1 n < < 6* 1 n 

C3^-lV3>/2 

C: x*e A < ClnC3>-l/'6>*x A 2+l Y2 V3+y 1 <x) 

C: deSoluety 1 =2*y/x*ln(3*yYx> , x, y> 

Type: Text 
Folder: 
Uariable: 
CEnter=0K 

wain-* 

choice__ 

ncdonald 

spit-al 

way el 
waoetxt 

hceD 

MAIN RAD HUTD DE 

MAIN RAD HUTD DE 


save copy as 


Type: Text 
Folder: mainj; 

Uariable: Idg 1 SI 


H 


■:: ;Eriler= 5 HUE; ::' CESC=CRNCEQ 


, e (2*Pl>*x 2/2+l/2>/3 I ans<i> 


[IE 12/30 


Im Texteditor wird diese Datei wieder geoffnet. Jede mit C: eingeleitete Zeile ist ein ausfuhrbares 
Kommando, das mit F4 "exekutiert" werden kann. Es empfiehlt sich, den Schirm laut Muster zu 
teilen. 

We open the file in the text editor. Each line showing a C: can be "executed" in the Home Screen 
by pressing F4. I recommend to split the screen according to the pictures. 


MODE 

T FI T F 2 T Fi 

1 ' 

|Page l|Page 2 1 Pag a 

2J 

■* Split Screen 

TOP-BOTTOM* 

Split 1 Hpp 

Horned 

Split 2 Hpp 

Text Editor* 

Humber of' Graphs . . 

l*. v 

Split. Screen Ratio 

2l 1+"'" 

Exact/ Approx 

AUTO* 

▼ Base 

DEC * 

(TEnter=SRUE"> 

CESCaCRNCELS, 

USE 4 - AND + TU OPEN CHOICES 




ft F 2 ^ "Tf^-Y FH Y FE 

■ [Command |U i ew |Execute|F i nd... | 


, 2-z- ln(3-z) - 




ln(2 • ln(3 • z) - 1) 


= ln(|x|)n 


C:d(z<x)*x, x 3 

C:Taz<2*z*ln<3*z;>-z;>,z;>=/Clx>: : ,x,c;> 
C:|ln<2*ln<3*z::'-1 V2=ln<abs(x)>+c I z=y^x 


'FiTTrii F2 T t f:: T t fh Y F - T 

— |Command |M i ew [ExecuteJF i nd... | 


2 ■ z • ln<3 • z) - 


;fr-m 


ln(2 • ln(5 • z) - 1) 


= ln(|x|) h 


Ciiitz'ixi'+x, x ; 

snach Trennung der Uari abler 
C: / < lz<2*z*ln<3*z)-z) , z)=J l Clzx,x,c.> 


Diese Frage stellen wir DERIVE 6. Der eingebaute 
"schrittweise Simplifier" verrat uns die angewen- 
deten - elementaren - Rechenregeln. So konnen 
wir von einem CAS nicht nur die Ergebnisse ra- 
scher erhalten, sondem gewinnen auch Kenntnisse 
in den Regeln - oder wir werden wieder an sie 
erinnert. 


Derive 6 gives the answer. Stepwise simplification 
shows us the secrets how to find the closed form 
of the integral. 

So the CAS does not only help us finding the re- 
sults very quick, but also helps us to gain new 
knowledge or reminds us on possibly forgotten 
knowledge. 


# 1 : 

2 • z- LN( 3 • z) - z 

If x>0, 

LN(x.z) * LN(x) + LN(z) 

1 

# 2 : 


dz 


z-( 2 .(LN( 3 ) + LN(z)) - 1 ) 
F(LNCa.x)) 


dz 


# 3 : SUBST 


dx SUBSTCJ F(x) dx, x, LN(a-x)) 


dz, Z, LN(z) 


2-Z + 2 • LN( 3 ) - 1 
LN(a + b-x) 


dx 4 

a + b-x b 

f LN( 2 • LN( 3 ) + 2-Z - 1 ) 

# 4 : SUBST , z, LN(z) 


SUBSTCF(x), x, a) * F(a) 

LN(2 • LN(3) + 2 • Lll(z) - 1) 
#5: 


# 6 : 


LN(2- LN(3-z) - 1) 
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Beispiel 3/Example 3 

(^3 x 2 v + 8 xy - y + 2^dx + ( x 3 + 4x 2 - x + — dy = 0 

V 


a) Zeige, dass hier eine "exakte DGL" vorliegt, d.h., dass die Integrabilitatsbedingung erfullt ist. 

8P dO 

P(x,y)dx + Q{x,y)dy = 0 ist exakt, wenn gilt: — = — . Show that the DE is exact. 

8y dx 

b) Lose die DGL auf traditionelle Weise. Erst allgemein, dann fur y(0) = 2. 

Solve in the traditional way, general and special fory(0)=2. 

c) Erzeuge fur die Losung dieses Typs von DGL ein Programm. 

Produce a program for solving DEs of this kind. 

d) Erfinde eine eigene exakte Differentialgleichung (die Winkelfunktionen enthalt). 

Find an own exact DE containing trig-functions. 

e) Suche die Losungen von a) und b) auf direktem Weg mit dem CAS. 

Find the solutions for a) and b) directly using CAS. 

f) Lose auch deine Aufgabe mit deSolve und mit dem Programm. 

Solve your problem (d) with deSolve and with the program. 

g) Zeichne das Richtungsfeld mit der speziellen Losung aus b). 

Plot the direction field including the special solutions from b). 


a) 




■-|j(3-x 2 y + 3x-y-y + 2)=-£;(x 3 + 4xA 

true 

... v-v+2 . y>=ri<x A 3+4x A 2-x+ y/2 . x>| 

MAIN RAD flUTD DE 1. 

b) Der allgemeine Ansatz fiir die Losung u(x,y) lautet: u(x,y) = j P(x,y)dx + /'(>’). Aus dem so 

gewonnenen u(x,y) wird nach Ableitung nach y ein Ausdruck fur r'(y) bestimmt, der dann nach y 
integriert zu r(y) ftihrt. 


f^lfllSralcSfJofflyp^Iolci^ UpI 1 

l?*a ImSralciTJoffl’erlp^Iofcle^ UpI 1 

1 qebraLaibafotherlprqm I olciean UpI 1 

■ + 4 ■ - x +-j + q(x , y) Done 

■/p(x , y)dx + r(y) + u(x, y) Done 

■ u(x , y) x 3 - y + 4 - x 2 ■ y + x (2 - y) + r(y) 

■-^Cu(x,y)) = q(x,y) 

-^(rCy)) + x 3 + 4x 2 -x=-^- + x 3 + 4x 2 -x 

■Jj^jdy + r(y) Done 

■ u(x , y) = c 

2 

x 3 -y + 4-x 2 -y + x-(2 - y) +-^p = c 

2 

■ ^4-x 2 ■ y + x (2 - y) + -^- = c | x = G and y = 2 

1 = G 

2 

■ x 3 y + 4 x 2 -y + x (2 - y) +-^p = 1 

2 

x 3 y + 4-x 2 y + x(2-y) +-^- = 1 

d(LL(x.v) . V> = Cf<X- v>l 

uXx.v>=cl 

...*3*v+4*x A 2*v+x»<2-v>+>j A 2/4=l 

MAIN RAD HU TO DE 6730 

MAIN RAD HUTD FUNE 11730 

MAIN RAD HUTD FUNE 11730 


Von den beiden Funktionsgraphen gilt natiirlich nur der obere! 


1 x i q |br alca Ic-Ioiharlpr^M 1 0 Ic 1 ealT UpI 1 

PiSfeK, A:a|Fa4Ap,-l -v-.!A t v R 1 

■ solue|x 3 y + 4 - x 2 • y + x (2 - y) + = 1 , yj 

y = 2(jx S + S-x 5 +14-x 4 -8x 3 + x 2 -2x^ 



■2(jx 6 + 8x 5 +14x 4 -8x 3 + x 2 -2x+l^ 

Done 

"\ , 

3 

-2*x+1 >+x#< x A 2+ 4»x-l > > 2 ( x>l 

MAIN RAD HUTD FUNE 11730 

MAIN RAD HUTD FUNE 


Only the upper graph is 
valid. 


c) Ein Programmvorschlag/Proposal for a program (following an idea of Gunter Redl). 


1 - 1 Icontr-o 1 ifxoluar If iind.. Jllcfdel 1 

(^bl^WATJo^plp^IoIaa^ UpI 1 

1 ■! Id = A, : CVU'-.cLh ’ ,J i 

:exakt(p_,q_7 
: F'r gn 
: Local u 

: I f d (. p_ , y ) *d <. q_ , x ) Then 
: Disp "nicht exakt!" 

: Goto end 

:EndIf 
:XCp_, x)>-Hj 

: u+XCq_-d<u, y) , y)=c+lsg_gen 
:Dis.p lsg_gen 
:Lbl end 
: EndPrgn 

"1, x ■ y • cost 2 ■ x) - cost y ) 

■ deSolue y = oaat2 . x:i x . 6l nC2->a ► 

4 1 2 X E 

4 ' X 3 k 

cos(2 ■ x) • y + 2 • x ■ sin(2 ■ x) ■ y + — ^ 4 ■ x * 

■ exakt|^3 ■ x 2 y + 8xy-y + 2,x 3 + 4-x 2 ->> 

Done 

2 

x 3 -y + 4 x 2 - y + x (2 - y) + = c 

2 v+8x*y-y+2 . x' 1 3+4x' 2-x+v/2>l 

MAIN RAD HUTD DE 

MAIN RAD HUTD DE 1S730 

MAIN RAD HUTD DE 17730 


p 14 


A Differential Equations Workshop for Handheld 


D-N-L#54 


dP dO 

d) Ein Vorschlag ist, mit mit Kern — = = s(x,y) zu beginnen. Dieser Term muss sowohl nach x 

dy dx 

als auch nach y integrierbar sein und ergibt dann die Terme P(x,y ) und Q(x,y), die noch mit 
jeweiligen Integrationskonstanten erganzt werden diirfen. 

s(x,y ) = x • cos(2x) + sin y 

P(x, v) = I* s(x, y)dy + konst - x ■ y ■ cos(2x) - cos v + x 1 
z -®- Q(x,y) - j s(x,y)dx + konst = ■ 

(x ■ y ■ cos(2x) - cos y + x 2 ) • c/x + 1 ' u ' Joy ^ JL> + + x • sin y-y 2 ■ dy = 0 

e) und f) Die CAS-Losungen: 


cos(2x) x-sin(2x) . 

— + - — - + x • sin y - y~ 

2 

cos(2x) x-sin(2x) . 2 

: + ] - + 1 — - + x-sm y-y 




iva aviAri:,,.; : | 


* J 1 x — 1 y t o 1 x 1 y yT^ZTx — t h • x — ■ J- x^ 

Done 

i exakt^x - y ■ cos(2 ■ x) - cos-(y) + x^ , cos ^ 

Done 

1 lsg.gen 

cos(2-x)-y . x-sin(2-x)-y , x 3 

4 + 2 + “ 3 — x ' c-osi;^ 


-2 Nx 6 + 8 x 5 + M-x^-S x 


3 + - 2 x + ► 

Done 


-I 3 ■ x^ • y + 8 ■ x • y - y + 2 1 

y T —1 0 ’*’► 


2-x J -y + 8-x^y-2x(y-2) +-^~ = G1 


„,y-y+2 )/<x 3+4x 2— x+y/2 > , x, y> 


y' =■ 


■( x ■ y • cost 2 • x) - cost y ") 


cos(2-x) . x-sin(2-x) 


oost2 ■ x) ■ y + 2 ■ x ■ sint2 ■ x) y 


2 

4 ■ x 3 


- 4 ■ x 


. *x A 3/3-4*x*cos <y >-4*y A 3/3=P2 


. >*y/2 + x A 3/3-x*cos<y>— y A 3/3=cl 


Im dritten Bild wird die Programmldsung der selbst erstellten DGL gezeigt. Sind die Ergebnisse 
mit deSolve() und dem Programm Equivalent? 

Are the deSolve()- and the program solutions equivalent? 

g) Das Richtungsfeld mit der speziellen Losung. 



y \\ \ / / / 

l \\\ / h 

/ i V — ""Cf 

\ \ \ \ \ \ \ \ \ \ 
!} \ IHUSU 
1 I \ W\ ws\ 

V s v-*'- 

\ \ \ 1 

\ 1 / -S.S V \ 1 

■| 1 f S s. V \ 1 1 

\ 1 //^\\' 1 s 

i / / ~ 

" / / 

■- \ i ////// y 

v i j i / / / / / / 

MAIN F:HD HtlTD DE 



Fiir „unsere“ DGL (d) enstehen recht barocke Losungskurven (ncurves = 4). Leider kann man am TI 
keine implizit gegebenen Relationen (Funktionen) zeichnen lassen - mit Umwegen iiber den R 3 ist 
dies recht muhsam moglich - siehe spater. Uber den Funktionseditor wird entweder das EULER- 
oder das RErNGE-KUTTA-Verfahren fur eine numerische Approximation eingesetzt. Ich 
vergleiche die Ausgabe am TI mit der Ausgabe von DERIVE. 

"My personal DE" from (d) results in very baroque solution curves. Unfortunately we cannot do 
implicit plots on the TI in an easy way - we can do it via R 3 representation as shown later. I com- 
pare theTI-plot with the respective Derive plot. 


Pfffes-jEffitrwaaiAS*..! l 

l?i2fc^Wlc,lR.<?«*liShlDK' u r2 SW 1 

i?i2fcs n wic,fe.<?^hShiDK' u r2 i 

jf’LDTf 

tO=Q. 

-(t • y 1 • cos(2 ■ t) - cos(y 1) + t 2 ) 

yl cos(2 ■ t) . t-sin(2-t) . , . y ,, k 

4 + 2 + i ■ sintylj y 

yii=" 
y?' = 
yi2= 
y3' = 

■-n 

yil= 

y / / — 

! / / ^ s s — 

I ; 

y 7 i / / 
01 Mil 
1 \ % \ \ 

1 \ \ M 

pflr. 

1 11/ rs /Sp-E'- 
i 1 1 1 / /tf - • V 



\ \ w i p l\ 

T \ M 1 pp — 

i i i i \}£- 

//it 

N - / 1 \ 1 1 

/ ) i 

^ ^ / t i 

— — / t i 


SELECT ONE 1ST-0RDER FUNCTION ONLY 

MAIN RAD HUTD DE 

MAIN RAD HUTD DE 
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Eine Verbesserung der Darstellung lasst sich durch Verkleinerung 
des tstep -Werts bei gleichzeitiger Verkleinerung des Bildaus- 
schnitts erreichen. Dass man auch iibertreiben kann, zeigt die Ab- 
bildung rechts oben. (Zooming in does not necessarily improve the 
graphic representation.) 


Im linken Bild sieht man die DAVi/KA’-Darstcl lung der implizit ge- 
gebenen Losungskurven mit den An fangsbed i ngungcn y( 0) = -1.8, 


Beispiel 4/Example 4 

2y" + 3/ - 27 v = 0; y(0) - -2, y'(0) - | 

a) Bestimme die graphische Losung am TI. Graphic solution on the Tl. 

b) Lose die DGL auf traditionelle Weise - CAS-unterstiitzt, aber ohne desolveQ. 

c) Vergleiche mit der reinen CAS-Losung. Compare the traditional with the CAS-solution. 

d) Interpretiere die Werte in der TI-Tabelle ([TABLE]). Discuss the [TABLE]-values. 

e) Ein Richtungsfeld ergibt sich iiber F9 Format Fields DIRFLD im [GRAPH]-Fenster. Versuche, auch die 
Losungskurve zu zeichnen und vergleiche mit den friiher gefundenen Graphen. Plot the direction 
field together with the solution curve from above. 


a) Fur die Eintragung im Funktionseditor muss ein kleiner Trick angewendet werden. Die vorliegen- 
de DGL muss zuerst durch geeignete Substitution in ein System von zwei DGL umgewandelt 
werden. The given DE must be transformed into a system of two DEs. 


Man setzt: yl = y, y2 = y', y3 = y" usw. Durch Differenzieren ergibt sich dann: 
yl' = y', y2' = y", y3' = y"' usw. 


Bei der vorliegenden Aufgabe bedeutet dies: 


„ 3 , 27 

y = — y + — y 
2 2 


Fiir die Zeichnung der Losungskurven miissen iiber F7 
die Axes nach der Abbildung eingestellt werden. 

Format Fields FLDOFF 



J-LDTJ 


^yl '=| 

f HUES ^ 


y i l = ' 
y2 1 = 
■-ii 2=1 
y3' = 

yi3= 

y4' = 

Ml 4= 

Fixes.. CUSTOM + 

X Axis t* 

V Axis yl+ 

<^Enter=SflUE) CESC=CflHCEL5 



vl 1 <0=v2 

USE t HND + TO OPEN CHOICES 


zAt|‘- |,ili rC'dnVi... 1 

JLDTS 

Plot is 
tQ=Q. 

^yl 1 =y2 
yil=-2 
y2 1 = ~ 3z2 
Mi2=lz2 
y3'= 
yi3= 
y4‘ = 

yi4=B 

m 2 + 27/2 yl 

vi4= 

MiilN 

RHD liUTO DE 


lOTb^Wic.fe.^liShlDfSwF/ 1 

Ml 



\ 1 

MAIN RHD HUtn DE 


b) Wir losen die charakteristische Gleichung, die zwei reelle getrennte Losungen liefert. Die allge- 
meine Losung wird als Summe von zwei Exponentialfunktionen mit unbekannten Koeffizienten 
angesetzt. Diese werden leicht aus einem linearen GLS berechnet, das aus den Anfangsbedingun- 
gen gewonnen wird. Mit geeigneten [WINDOW]-Werten lasst sich die Losungskurve schon darstel- 
len. Solve the characteristical equation and find the general solution as usual. 
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ft 1 |r 1 qebra Ic^L: lother IpTfn I olcieaif UpI 1 



( 9 ■ @2 



A = 3 or A = _ 9/2 

■ i?l -e 3x + @2 ■& " ^2 x + Done 

■ b4CG)= -2 (£1 +62 = -2 

■-^ b4<x M = 1 '*2|« = 0 3-ei - u -2 

61 = - 17/15 and 62 = - 13/15 
■ 61 ■ <? 3 ’ x + 62 e " 9/2 ' * | 61 = - 17/15 and <!► 

-17 -# 3 ' x 13-e 2 

15 15 



c«vi<x>.x>lx=0 

xVi5-i3*e' i <-9»x/2Vi5+vi<x>l 

MAIN F:HD fiUTD FUNC 9.' 9 

MAIN F:HD fiUTD FUNC 9Z30 

MAIN FiAD HUTU FUNC 


c) Die reine CAS-Losung ergibt sich sofort in allgemeiner Form. 
Zur speziellen Losung muss man vorgehen wie in b) gezeigt. 

We obtain the pure CAS-solution in one single step. 


~y F3T Y FHt Y F5 Y FfiT Y ) 
a|ualc|Qther|PrgnIQ|clean Up] | 


■-^CyiCx))|x = 0 
■ deSolve(2 y 1 ' + 


1 b4(x) 



1/2 

- 27 ■ y = 0 , x , y) 

-9x 

= i!2-f 3 x + ei f 

2 


"9-x 

®t‘f 3 ‘ x + e2-s- 

2 


b4<x>l 


d) Die Tabelle wird zuerst eingerichtet, dann kann sie betrachtet werden. 
Unter yl und y2 werden die Werte fur y{x) und y'(x) ausgewiesen. 


i^u..rA. it, ; . l 


^feKupMUlH,:;: ; ,Ti:Tp«krT-F,,f 1 




f TABLE SETUP 



t 

y: 

y2 




t 

J=L 

y2 




tbl Start EH 1 

. 1000000 

-2.Gb998 

-2. 1G839 


-. 30GGGG 

-3.70042 

13.24332 


atbl 1.1 1 

. 2G0GG0L 

-2.39201 

-4.58792 


- . 200000 

-2.70645 

7.556784 


Independent.... FIUTO + 

. 3000000 

-2. 98884 

-7.27283 


-. 10GG00 

-2.18165 

3.549339 


. 4G 00000 

-3.83438 

-10.4726 


0.000000 

-2.00G00 

. 5000000 


. 5 G 0G GOG 

-5.05513 

-14.5151 


. 1G00GG0 

-2.06998 

-2. 1G839 


CEnter=SflUE!> CESC=CflNCEL} 

. 6000006 

-6.73186 

-19.7853 


. 2000000 

-2.39201 

-4.58792 


>uuuuuui -y . uuyubi -zt. . zbbzi i 


. 7G 00 GOG 

-9.G08O5 

-26.7652 


. 3G00G0G 

-2.96864 

-7.27283 


t 

=□. 


a 

II 

JJ 


1 

II 

TYPE + £ENTER1=DK AND IESC1=CANCEL 


MAIN FiAD HUTU DE 


MAIN FiAD HUTU DE 


Setup the table and inspect the values of y(x) and y'(x) as yl and y2. 


Wenn man auch die Werte der zweiten Ableitung sehen will, dann muss man wieder in den 
FUNCTION-Modus zuriickwechseln. Der Vergleich zeigt den „Fehler“ bei der hier angewendeten 
EULER-Methode. Wenn man im DE-Modus auf die RUNGE-KUTTA-Methode umstellt, dann 
fallt der Fehler deutlich geringer aus. 

For finding the values of the 2nd derivative one has to switch back to FUNCTION-Mode. Compa- 
rison shows the "error" applying EULER-method. Applying RUNGE-KUTTA results in smaller er- 
rors. 


IGilAT eIAi ylT.CUhA 1 




v^, r :u iV- , : 1 

JLDTJ 


X 

y: 

y2 

y3 



t 

y: 

y2 



, -17 f 3 * 13 ■ f 2 

IDHDIDIDIDIDIO 

-2.00000 

. 5000000 

-27.7500 

unnlnlnmlnin 

-2.00000 

. 5000000 


. 1G0GGG0 

-2.G8245 

-2. 1G277 

-24.9589 

. 1G0GG0G 

-2.G8236 

-2. 1G296 


-y2^yl(xa 

d 2 

. 2000000 

-2.41743 

-4.6G958 

-25.7289 

. 2000000 

-2.41717 

-4.6G984 


. 3G0GGG0 

-3. G 1223 

-7.35161 

-29.6376 

. 3000000 

-3.01162 

-7.35142 


. 4000000 

-3.906G6 

-10.6437 

-36.7662 

. 4G0OG00 

-3.90469 

-10.6419 


. 5G 0G G00 

-5. 17G59 

-14.8267 

-47.5630 

. 500GG00 

-5. 168G4 

-14.8214 


dx 

y4=B 

. 6000000 

-6.91451 

-20.3067 

-62.8859 

. 6000000 

-6.91G11 

-20.2956 


. 7G00G00 

-9.29213 

-27.5979 

-84.0470 

. 7G00G00 

-9.28474 

-27.5774 


v4<x>= 


x=0. 


O 

II 

MAIN FiAD HUTU FUNC 


MAIN FiAD HUTU FUNC 


MAIN FiAD HUTU DE 


e) Im Richtungsfeld wird der Zusammenhang zwischen yl (= Losungskurve) und y2 (= erste Ablei- 
tung) dargestellt. Stellt man die Achsen auf TIME um und aktiviert yl' und y2', dann werden Funk- 
tion und erste Ableitung gezeichnet. (Flier wurde die Funktion mit dem Style Thick ausgezeichnet. 
AnschlieBend habe ich die Graphen der exakten Losung und deren erste Ableitung fiber die nume- 
rische Approximation gelegt (RUNGE-KUTTA). 


i'FlTTTlY F2TY F3 FH Y FST Y FhT YF? aY FB T 

| -■ ■ Z oo Di | T ra ce F: e gr aph |M at h | Dr a w W V 1 1 C | 

i?i2fcs , n wic,fe.<?«*hShiDK' u r, ? sw i 

I^IzSTT^cdReSaPhlnLYhlD^J?^!^ 1 



y 


S 

» 

n ) } ) ) ) 1 1 ) / 
!///////// &*<--- 
////// 

^ \ \ \ \ 
N \ \ 
K V \ 
Y \ 

222x4 

< 

\ • 

r 


MAIN FiAD HUTU DE 

MAIN FiAD HUTU DE 

MAIN FiAD HUTU DE 


Wird die EULER-Methode eingesetzt, kann man leichte Abweichungen auch optisch erkennen. 
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Beispiel 5/Example 5 

y' - 3 / + 10 v = x 2 - 2; y(0) = 1, y'(0) = -1 

y( 1) = 1, y(l) = 2 

a) Bestimme die graphische Losung am TE Graphic solution on the Tl. 

b) Lose die DGL auf traditionelle Weise - CAS-unterstiitzt, aber ohne deSolve(). 

c) Fiihre die "Probe" durch. First without and then check with deSolve(). 

d) Vergleiche mit der reinen CAS-Losung. Compare with the pure CAS solution. 


e) Die beiden Losungskurven werden leicht gefunden. Fur die zweite sind die Anfangswerte zu ver- 
andem. Diese konnen auch interaktiv liber F8 im Grafikfenster eingestellt werden. 


l^lzSlEdid^lH P RIsaiJ R S S J 1 

fFiTfUDT FST T F3 Y FH Y F5T Y Ffi’r YF? & 

t t — Zoom T race Regraph Nath Draw ▼ V 

F0 1 

ID 


j-ldtj 
tu=ai 
^yl 1 =y2 
yil=l 

^y2 1 =3 - y2 - 10 yl + t 2 - 2 
y i 2= " 1 
y3' = 
yi3= 

y4' = 

yi4= 

td=o. 

y 






y 

. Ay 



r LAA 

J 



L 


~ l\ 1 

t 

1 

CHECK HUES AND INITIAL CONDITIONS 

MAIN F: AD AUTO DE 

MAIN HAD AUTO DE 


Besonders interessant wird 
hier das DIRFLD fur Funkti- 
on und Ableitung. DIRFLD 
for yl and y2. 


sp 

■nax=5 . 


^nax=b. 
^step=. 01 
L.plot.=0. 

<n i n= "50 . 
criax=2GG . 
<scl=50. 
jm i n= - 60 . 
jr-iax=32G . 
jsc 1 =5Q . 
icuryes=0. 

isiep^ l . 
i~ ii~P=.= 'T I . 


DE 



■y2f / / 
/ / / / 
/ / / / 
/ / / /j 
/ / /// 
/ // / 
/ /// / 

S S S ^ N N \ \ 

^ ^ ----- \ \ \ \ 

\ \ \ \ \\ \ \ 
\ \ \ \ \ s u\n 

> \ \ \ \ l c \ \ \ C 

\ \ M M M 1 1 1 III 1 

\ v-— 

y ///////// \ \ \ \ h \ 

MAIN RAD M U T □ DE 


b) Die traditionelle Losungsmethode beginnt wieder mit der charakteristischen Gleichung, die hier 
zwei konjugiert komplexe Losungen ergibt. Der Realteil wird fur die Exponentialfunktion und 
der Imaginarteil fur den trigonometrischen Losungsansatz verwendet. 


If 1 2“ Im 1 qebralcalc-Ioiher Iprcfn I die lealT UpI 1 

PiSfei ^EraTcH’cToiSrbrSnlofcl^n UrI 1 

I t 1 ^2 Ifl 1 qebraka Lclot her Iprqro I ole 1 ean UrI 1 

■ solye( A 2 - 3 ■ A + IQ = Q , a) fa 1 se 

■ cSolwe( A 2 - 3- A+ 10 = 0, a) 

x = Z/2 + ^- i or- > = 3r2-J|i- i 

■ a (ti sln (M | , co ,|M . || . ,► 

■ 4 | e 1 si n[ x J + S 2 cos^^ ■ x J J a y 0 ( x) 

Done 

■ bO + bl ■ x + b2 ■ x 2 + st(x) Done 

■ -^3— [ stlxl - 3- -Ac stt X)) + 1 0 ■ stfcxj = x 2 -► 

ax 

10b2x 2 +(10bl-6b2)x+10b0-3blF 

■#{ 10 -b 2 =l and 10 ■ bl - 6 - b 2 = 0 and 10 -!> 
bO = --ift and bl =3x50 and b2 = 1x10 

■ b 0 + bl ■ x + t .2 ■ x 2 I b 0 = - and bl = 3X5> 

Done 

■ y 0 ( 0 ) + stCO) = 1 C 2 -|H’-I 

,*x> + e2*cos(J(31>/2»x)>+vD<x>i 

*1.2 > *x+ lG*bD-3*M +2 *b2 = x*2 -21 


MAIN FiAD AUTO FUNC ifi^O 

MAIN FiAD AUTO FUNC ifi^O 

MAIN FiAD AUTO FUNC lfir'I-Mj 


Dabei bleiben zwei Parameter. Fur das Storglied x 2 - 2 wird ein unbestimmtes quadratischen 
Polynom st(x) angesetzt, dessen Koeffizienten nach Einsetzen von st(x) in die DGL durch 
Koeffizientenvergleich gefunden werden (bO, bl und b2). Durch Einsetzen der Anfangsbedingun- 
gen erhalt man schlieBlich die noch offenen Parameter. The traditional way supported by the Tl. 


d 

dx 


'A-(y0Cx) + st(xJ)= - 1 | k = 0 


e l. E2 -W = 


ei'JST 3 ■ @2 
2 2 
si m , 


1 and 

ei = 


2 

2865 J5l 


5-02 

2 


+ 3^ 


15500 Jriu ^ 5 M 

... 1 >/2+ 3 *02/2+3/50= ~1 . <01 . P2> >1 

MAIN F:HD hUTD FUNC lb.- 1 SC 
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Das ist die Losung fur das erste Paar von Anfangsbedingungen - in „Schonschrift“. 

601e T cos^* 2863 V3le T sin 2 . 101 

. . 9 9 x 3 x 101 

v(x) = — + — + 

500 15500 10 50 500 

The solution for the first pair of initial conditions. In c) we check the result, 
c) Die Probe bestatigt unser Ergebnis. 


" 1 < ^ + ~T5=L 


Done 


2(yl<>0) " 3~^(ylCx)) + lG-ylCx) 


■{yl(x> -^(yHx))} | x = 0 


d) Wir losen die Aufgabe mit desolve() und verwenden die allgemeine Losung zur Berechnung der 
speziellen Losung fur das zweite Paar von Bedingungen. 


upT i upT i m \ \ s u 

H X H I l, jZ J K 2 . } ^ ~~ 


'(yiCx) -^(ylCx))} |x = G Cl -1J 

1 deSolyeCy 1 ' - 3-y 1 + 10 - y = x 2 - 2 , x , y) 

3 ■ x , 3 ■ x 

y-Gl-f 2 cos 


1=0! .'Jfifl ,3X2 + e2 . sin (J|L> 




82 t 2 - siril — ^ 


fJUl 




'■i e2 = 


177 cos[J|ti] f ■ 5/2 - l3l 521 -si^, 
15500 * 


... > *x/2 > ± < 5 B*x^2*30*x-1 Q1 > /5 OQ ...21*sin<f <31>/2>*e*<-3/2>/SQa ... 21*sin <J<31 V2 >*e*-< ~3/2 >/S0B 

HAIM F:HD HUTD FUNC 17.-' SO MAIN F:HD HUTD FUNC Sii.'jii MAIN F:A[i HU TP FUNC l/ZO 


Die komplette zweite spezielle Losung kann mit einem mit DERIVE produzierten Ausdruck ver- 
glichen werden: Compare the 2nd special solution with the respective Derive result. 


DS01VE2_IV(-3 , 10, x - 2, x, 1, 1, 2) 


3-x/2 





f /31-x 5 


cos 



V 

l 2 J 

\ 


- 3/2 
521- e -COS 




- 3/2 

177-./3l-e -SIM 


500 


15500 


' ,/3l ' 

A 

„ 2 . 


> 


f - 3/2 

177-./3l-e 


,/31-x 


Sill 


H-t-) 


15500 


- 3/2 

521-e -Sill 


500 


.31 ’t 


2 J 


50- x + 30- x - 101 


500 


Beispiel 6/Example 6 

Eine Anwendung: Schwingungsgleichung Vibration Equation 

Die Leder eines PKW wird durch das Gewicht des Lahrzeugs (ca. 3500 N pro Rad) verformt. Die Le- 
der wird um weitere 5 cm verformt und dann sich selbst iiberlassen. Die Schwingungsgleichung s(t) ist 
aufzustellen und der Graph (Weg-Zeit-Diagramm) zu erzeugen. 

Dabei ist m die Masse (m = 350 kg), k die Lederkonstante (z.B. 50000 kg/s 2 ) und / der Dampfungsfak- 
tor. Verschiedene Werte fur / fuhren zu verschiedenen Lallen. Untersuche die folgenden Lalle und 
interpretiere die entsprechenden Graphen: 

a) Dampfungsfaktor / = 0 dvarious damping factors 

b) Dampfungsfaktor / = 1000 kg/s 

c) Dampfungsfaktor / = 10000 kg/s 

d) Welcher Dampfungsfaktor ergibt eine gerade ausreichende Dampfung? 

Which damping factor is necessary for sufficient damping? 
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Alle Aufgaben werden mit grdBtmoglicher CAS-Unterstiitzung gelost. 

This example was presented by Gunter Redl in a T3-ACDCA-Seminar. Many thanks, Gunter. 


a) Es ergibt sich das Bild einer ungedampften Schwingung. undamped oscillation 


l?«=feu3Sr : ,T.|-A.icl,:. i r • 1 

I t 1 ^2 Ifl 1 qebralca fclot her Iprqro I O lo 1 ean UpI 1 

l?^bS,Wlc.lR.<5?«*klShlDS’ u r. ? 1 

■ deSolveCn • s 1 1 + k ■ s- = G , x , s) 

f 10-470 -x) . „ . [' 10-470-x 1 

s- = G1 cos| j 1 + 122 sinl ^ 1 

_J ( 10-470-xl . . f 10-470 x1 . , , 

■ T cos[ =; 1 + G2-sin[ j 1 + s(>0 

Done 

■ S-CO) = .Q5 li^mBUnlB 

■<.-cos[ 1U 'V IJH J + e2-sin[ iU ' J ; l:JX J^Cx) 

Done 

■ sCO5=.05 ei = .Q5G0 

-&«M)-0|*-0 !'• - 0 

■ . 05 - cos^ 10 X j ■* y 1 ( x) Done 

\ / 

\ A A A A A 

l 1 

v v v V/ V v 

0 . 05 *c o s <1 0*J < 70 > »x/7 > -» vl < x > 

H . H5*ooi; < 1 0*1 C 70 >*x • 7>->v1 < x >1 

MAIN F:HD HUTU FUNC 3Z E 

MAIN F:HD HUTU FUNC 6Z30 

MAIN ft AD HUTU FUNC 


b) Hier sehen wir eine zu geringe Dampfung - eine gedampfte Schwingung. Die Funktion s(t) ge- 
meinsam mit ± 0,05e‘ (10/7)Y . too less damping 


1^1 -.At-A-aLf ’ . 1 

1= |r i qebr a fcaT:- [other |pr-gn I olc-lean UpI 1 

I^I^TEeU.Arhh^'.^TV,. 1 

■ deSolueCn ■ s 1 1 + 1000 - s' + k ■ s = 0 , x , s) 

-10-x . _ .. -10-: 

s = 07-a 7 f 7 > 

-10 X 

u 7 7 J 

■ solwej 10 ®S‘J59 _ is_®7 =0,es]|07 = .05 

es = . 0060 

■ s<jO | G7 = .05 and @8 = . G06D192926542SS5 

-lG-x 

< :,J f • .9060 . 7 - si n(‘° f *] 




. 4-469 -xl . „„ 7 . f 10 -459 x1 . , , 

a ~f J + G8-# -sinl ^ 1 ■* s(sO 

WRIT* 




ans <1 >Ay2<x> 

ans <1 >Ay2<x>| 

MAIN FiAD HUTU FUNC b.'lH 

MAIN FiAD HUTU FUNC lH.-'SD 

MAIN FiAD HUTU FUNC 


c) Mehr als ausreichende Dampfung - Kriechfall, aperiodische Bewegung. too much damping 


I^^SIr 1 qebralcalclotherlprqp I olciean UpI 1 

|^|^ AIT T.A-hA 1 


■ deSolveCrv s 1 1 + 1G0GG -s' + k • s = G , x , s) 

( 10 -430 100 1 ( -10 430 ^ 

s = 010- ft 7 7 J + 09 A 7 

1' 10 430 100 1 I" -10 - 430 

■SIS -ft 7 7 j 

i 10-430 100 1 r -10-430 100 1 

jldt; 

y2=.05-+ 7 -p= 1J E V % -■ 90.-0: -)?> 

y3=.05-e " 10 '' 7 ‘ x 
y4=-y3(x) 

f 10-430 100 1 

7 7 x ^ 

^5=. 070643546458767 ■ * ( } - 

y6= 
y7= 





ans Cl >Ay5<x> 

v5<x)= <-10*J<30>/7— 100/7>*x> 

MAIN FiAD HUTU FUNC £2/30 

MAIN FiAD HUTU FUNC 

MAIN FiAD HUTU FUNC 


d) Eine gerade ausreichende Dampfung ergibt sich, wenn der Dampfungsfaktor 1 so gewahlt wird, 
dass die charakteristische Gleichung der DGL eine reelle Doppellosung hat. 


|T?2|fli^,- a Tc^T0^r|pp^I0lci^ UpI 1 

1= 4 ^+ |r 1 qebr alcaT:- [other Ip, - qn I ole lean UpI 1 

I^ES^Tr^ &e Ifie^ a phIwS,ra^/I:: 1 

■ soloefn • X 1 - + 1 ■ > + k = U, >,) 

J 1 2 - 70000000 - 1 -[■! 1 2 - 7000> 

A 7GG Cl ■" 70G 

■J 70000000 + 1 1000-470 

■ deSoluetn ■ s 1 1 + 1 ■ s 1 + k ■ s = 0 , x , s) 

-10-470 X 

s = tau ■* + eiaj-ff 7 

■5(01=. 05 012 =.0500 

■ -£-(s<x>) = 0 | x = 0 and 012 = . 05 

011 - .5976 = 0 

■ sCx) | 012 = .05 and 011 = .5976143046672 

-10 -470- X 

( . 5976 ■ x + . 0500) ■ f 7 

l 

V. 

deSolue <m*s 1 1 +l*s 1 +k*s=0, x, s> 

ans <1 >Ay6<x>| 

MAIN FiAD HUTU FUNC £5f , 3A 

MAIN FiAD HUTU FUNC WStt 

MAIN FiAD HUTU FUNC 


Choose damping factor I so that the characteristical equation must have a real double solution. 


On the next pages you can find problems from DERIVE Newsletter #2 and #3 ( 1991). Then 
we demonstrated the use of the MATH-file ODE1 .MTH. Some of the then useful auxiliary 
functions are missed in recent ODE1.MTH, because programming made it possible to inte- 
grate them into one procedure. 

I recommend to compare with the updated contributions given in DNL#2 and DNL#3. Josef 



Examples from DNL#3 worked out on the TI-92+/V 200. Compare with the DERIVE procedures! 

Ex 1 : (1 + y 2 ) dx + xydv - 0 

Give the general solution and give the special solution containing P(3|2). 

Sketch the direction field, some integral curves and the graph of the special solution 
from above. 

(1 + y 2 )dx + x ydy = 0 — > l + y 2 +xyy' = 0 



Ca 1 c- Other Prgn 1 0 


' soluej y = — 2 ~ 1 ? yj 


DrauFunc J<45-x A 2>/x 


I Miii h fihd hutd fun-: r j 



xc:3.G000 


IHlilM FiliP liUTD 


c: 2. 6000 


fun-: 


C — x 

General solution is y 2 = — . We plot the direction field together with the integral curve passing 

the given point. 



ncurves from the [WINDOW]-settings does not work - (tO is temporarily set at the middle of the screen 
and initial conditions are distributed evenly along the v-axis). As there are no defined values on the y- 
axis ncurves must fail! ! We can add solution curves via DrawFunc in the [GRAPH]-window or we change 
xmax and xmin to shift the y-axis out of the center (ncurves = 10) 





















„ . - , , , , , .. -COT (x) -1/SIN (x) 

Derive version 2 showed the solution: y = 2 



Plotting both functions yl and y2 on the same axes shows the identity of the solutions. Then we see 
the direction field together with one branch of the special solution. 
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Ex 5 : The general equation expressing the relation between electromotive force U and current / in a 

circuit containing the resistance R and the inductance L is: 

U-L -I' + R-I 

Solve this equation first for a constant U = Uq. 

7(1=0) = 0 

Then for U 0 = 20V, R = 5 Q, L = 0. 1 Henry. 

Solve this equation, when U= U 0 sin(cot)! (eo = 40) 
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Ex 6: The DE y' — x + y defines a nice direction field. Let the Voyage 200 plot this field together 

with a family of integral curves. 



SELECT DNE lST-URDER FUNCTION ONLY 



Ex 7: Find the general solutions of the given DEs and verify them. Find also the special solutions 

containing the given points. 

a) xy' -4y + 2x 2 +4 = 0; P(\ | -3) 

b) y'cosx + ysinx = tanx; P(f|l) 

c) (1- x 2 ) y' = 1 + xy; E’(2|2) 


■ debolut-l..:::: • y 1 - 4 • y + 2 ■ + 4 = 0 , x , yj 

y = @4-x 4 + x 2 + 1 

■ y = @4 • x 4 + x 2 + 1 | x = 1 and y = "3 

-3 = E4 + 2 

■ -5x 4 + x 2 + 1 + y "5-x 4 + x 2 + 1 

■ x ■ -rr-( y) - 4 • y + 2 • x 2 + 4 = 0 true 

x*rf<v.x>- 4v+2x A 2+4=0 

MAIN FiliD E W ACT DE 


^ | H ! g e b r a [c a lclo the r]p r g n I □ [c 1 e an U p [ , 
■ deSolueCy 1 • cos(x) + y sin(x) = tan(x) , x , > 


2 • @5 ■ ( cqs( x)) ^ + 1 
2 • cos(x) 



(26 = ln(j3 + 2) + 2-J3 


| (26 = ln(j3 + 2 > 
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and the solutions: 


upT i 

■ daSolue^y 1 = ^ , X , yj 


h X • y - -=^- = G 1 


(2 ,2 

1 solve! ^ 2 ~ + x ■ y - -H>- = 01 | x = 1 and y = 2 ■? 


solve <ans <1 > I x=l and v=2.Pl> 

MiilH FiiiD liUTD FllHi 2/30 




FrgmlQlCleai 


iTuX~1 


arid y(lx2) = - 3x2h 


(f 2 y + 3 - i) = -[lnCx) 


> x _ ,2-y . 


... A <v/x>> and v<l/2>=-3/2 .x.v>l 

MiilH FiiiD liUTD PE 1 ? 0 




Clean Up 


a 


<? “ 3 (f 2 y + 3 - l) = -(lnCxJ + f- x - f 2 y + :* 

■ iolue(f “ 3 -(f 2 y + 3 - l) = -( lnCxJ + f x - e* 
= x(ln( -(a 3 ' lnCx) + a 3 ' ln(2) - l)) - 3) 


...x>*e^<v/x>-e- x <2*v>*ln<2>>.v> 

MHIH RHP liUTD PE £/P0 


Ex 10: Compare the integral curves of 

, v 2 -x 2 , , 2x v . . , . 

y = - and y =— — — intersecting in P(2 \-\ ). 

2xy y —x 

Plot both curves together with the tangents in P. 

I tried to superimpose two 
implicit plots using the 3D 
mode. On the TI 1 cannot plot 
the tangents. 


fF 17 *irii F 2^ T F3^ T FH t j FE T Fb^ 

| ^ 4 — |h l gebr-a |Ca 1 c lother |Prgm 1 0 |C 1 ean 

UpT 



— j— + 4y + 11^3 = 


AiFlj 

■ expand] 


+ 4 y + 11X3= -(x 2 - 

4)- 

h\ 



y 3 + 12-y + 11 = 12- 

d-3 

2 

■ x ■ y 

■[ "y 3 + 12-y + 

11 = 12-y - 3 - x 2 ■ y) 

12- 

y 



11 - y 3 

= -3 

■ x • y 

an s < 1 > 

-12vl 



MiilH 

FiliP liUTD FUNC 1 B 3 0 



^l^Jci^b&Wnlolci.^ UpI 1 

■ deSolve^y 

° V x -y ^ y C25= -l.x.yj 

? -X-C2-X-5) 


y 2 

■ deSo 1 vej y 

, y ,, and y<2) "1 , x, y] 

y - x‘ J 

-hj— + 4-y + 11x3 = -(x 2 - 4)-y 

... /<v A 2— x 

A 2> and v<2>=-i.x.v> 

MiilH 

FiliP liUTD FUNC ifi/30 






hr 

MAIN FiliP flUTD 3D 




Ex 1 1 : Find the curves which intersect the circles x 2 + y 2 = r 2 at angles of rc/4. 

Plot some circles together with a family of curves representing the general solution of the 
underlying differential equation. 


v — X 

According to [2] the problem leads to the homogeneous DE y' = 1 — — . 

y + x 



f^Im^bralcKofaffirTr^Iolc-l^n UpI 1 


, ln 

x 2 + y 2 l , - t X yil 

2 +2-tan-i[ x j > 

- 



.(ooiCS)) 2 ! J ► 

2 


n.2»ji>-ji>/2 = lnfcos<0>*r> + Pi 

Mlilh 


FiliP liUTD PE 2/3D 


P^ImyTbnalcffoyirenlpP^Iolci^ UpI 1 

■ deSol'.,'e| 

f 

C/ 2 )— w) 


bins Cl > 1 

x=r*cos<0) and v=r*sin„, 

MAIN 

FiliP liUTD PE i/30 




■m 


cm 


= lntx) + Intel) 

+ 2-e 


= ln(x) + Intel) 


2>+2»0>/2=ln<x>+ln<ei> 


pe 



p 
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'-'pi 1 




tta 


= InCx) + Intel) 

+ 2-e 


= InCx) + lnCGl) 


2>+2»0>/2=ln<x>+ln<Pl> 


DE 3/30 


gEgEKjil upT 1 


l t r 



[- 

, ln 

x 2 + y 2 1 

+ 2 ' 9 J 

x 2 . 


2 


- = ln(!x) + Intel ) 


= ln(x) + ln(0 


-ln(x 2 + y 2 ) 


ln(x)-0 = ln(x)+ Intel:) 


EWgBBHjBE14 HI Pgi] 


Pfflm^bna 

FI* j FH* T FE T Ffi^ j l 

Calc|Other|PrgnIO|Clean Up| | 

■ expand] | — — 
-ln(x 2 + y 2 ' 
■( - InCx 2 + y 2 

- JS - £ - + InCx) - & = InCx) + InCGl) 

(x 2 + y 2 ) , 

^ + InCx) “ 0 = ln(x) + * 

+ 2- InCx) - 2-0 = 2- InCx) + 2 ■ > 

) + 2 ln(x)- 2-0 = 2- ln(x) + 2> 

-ln(x 2 + y 2 ) - 2-0 = 2- ln(61) 

ans <1 >-2*ln <x>l 

MHIN FiiiD liUTD PE B / 3 P 




-ln(x 2 + y 2 ) - 2 S = 2 -Intel) 

■ - InCx 2 + y 2 ) - 2 0 = 2- Intel) | x = n ■ cosf:®)> 

- InCn 2 ) - 2 6 = 2 Intel) 

■ so 1 ueC - 1 nC n 2 ) - 2-0 = 2- ln(61) , r) 

- f ' e f -® . 
n = — xi — and El > 0 on n = ' and El 


Isolve (a.ns <1 > ,!»>! 


MiilH F;HD HUTD PE B/30 


We find a family of logarithmic spirals intersecting the cir- 
cles. As you can see we switch to polar coordinates. Perform- 
ing this task requires a lot of mathematical and technological 
competence. 




T^|2S*:T t r Inli rCial) 1 

;,..i ) 

J-LDTS 1 

| 61 = {.05 .1 .15 

l»2=-|j-^ | El = 0.05 .1 .15 

.2 .3 . ^ 

.2 .3 .► 

IE^r|ei = C.E5 .1 .is 

2 .3 .41 

n4= I 

t-5= 

ir3<e>=i/eiiei={.G5. .j 

.-15..2....I 

1 MAIN FiliP liUTD F'DL 1 



Ex 12: A point moves on a curve in the x-v-plane in such a way that the angle formed by the tangent 

of the curve with the x-axis is three times the angle between the radius vector and the x-axis. 
What is the Cartesian equation of the family of curves satisfying this condition. 

Plot a family of solution curves, especially this one which contains P(-3|-2). 


fFi-THIOT FZ* T F3 j FH T FS* T Ffi’r YF7 
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MHIN 




:HP 

liUTD 




PE 






PiSfelSrartS'obfflrWS.IorCleSS’ UpI 1 

■ deSolve 

-2-(x 2 h 

■ ! 

Jx-J" 

, 3-x 2 "y-y 3 1 

y 'y- 3 ,,2’*' s J 

Jx-Jy 

y 2 ) . h 

= @ 1 x = r ■ cosC^) and y = r ■ * 

...Gi 1 x=r»cos(9> and v=r*sin(0>l 

MHIN 

RHP HUTD PE 3/30 



The first - graphical - approach is not really satisfying. Again application of polar coordinates proves 
to be a successful means to obtain better results. More explanations can be found in DNL#3. 


i„x, A i, i^| Hl y.T -A : ki.'lfnp I :• 2 J 

i * • .1 u ^ ... \ - A 


-2-r £ , 

JsinCG) • r •JgosCCit 


-2-|r| 


■JsinCS) • r • cosCG) -r 


= ei 


■2-|r| 


J sinCtO ■ cosCy) 


"2 ■ J2 ■ |r| 

1 solve ■ . == ■ = iif 1 , 


"I, -J s i nC y) - cost; y> ) 


-2-J2-11-1 

JsinC2-0) ' 


-2 J2 Irl 
Jsin(;2§:) 


.!•£. ! ni.rx' ■ GosCgj 


I = : x2 O' ' 


, , "2 ■ J2 ■ |r| 

1 so 1 ve - . . . ==- = e 1 , r 


I I = ix2 O' 
j 2 -=.it-n :2 e:i -ei 
4 


and Jsin(2-0)-El S 0 t> 


!■ liUTD PE j /30 


SQlue<ans<l>.r> 


s q 1 v e < an s < 1 > . r >1 
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^|zL, :T.r I'l VTTLL ] 

(^^Iz^lT^oefe^Phk^lp^l^ fW 1 

Ifi 1 qit.i-.alc ale fotherlp™ I ofc: 1 ean Up! 1 

JLDTS 1 

<rl ,J2-ainC2-e)-ei. | | B1 , I2 4 6 s 10> 

r2= 

r3= 

r4= 

r5= 

r6= 

r?= 

rS= 



■ 2 L ^ = (21 1 x = -3 and y = -2 

•lx -ly 

-2(x 2 + y 2 ) 13j6 

-Jx fy 3 



ri <e>=... .4.6.8.iG.12.14.16.18> 



MlilN FiliP liUTD FDL 

MlilN FiiiD liUTD FDL 

MlilN FiliP liUTD FDL 7/SO 


The solution curves are lemniscates. 



Ex 13: Show that the given DEs are of exact form. Give their general solutions and to each one the 
special solution if there is a point given: 

a) lx 1 - y 2 + v - y' ■ (2xy - x- 4y) = 0; P(l|-1.5) 

x t 

b) 2x + e x In y + = 0 

y 

c) 2xy - (2x 2 + y)y' = 0; Q{- 2 1 3) 



m - ■ - I ;1 - LI 

■ ^2-x-y - x + 4-y) = 0 and y(l)= _ 3^2,x,y) 
y + 33/4 = 2 3 * + x(y 2 - y) - y 2 + y + ^ 



... < an s < 1 > >-rixht Cans <1 > >=Q>*12 

MHjH FiiiD liUTD F PL 1.-' H 




^I'lPrSiIofciiwin ;Jr-l 1 

■# 2 '^ y 5-^C- 2 

x 2 - y) EBB 

■ deSo 1 ue( 2 ■ x ■ y - y 1 

(2 x 2 + y) = G, x, y) 


X 2 1 

i T + y = B2 
y a 

■ deSol ','■=•( 2 ■ x ■ y - y 1 

(2 ■ x 2 + y) = G and y( -► 

ans<l>»9y A 2 

MlilN FiliP liUTD 

F PL H/ ? 


Pi=|ni<aEr.lcK'c WKrWS’n.Iolci UpI 1 

■ deSolve(2 • x 

y - y 

■(2 "jc 2 + y] = 0 , x , y) 
x 2 i 

^2 + i = ® 2 

y * 

■ deSolve(2 • x 

y - y 

■( 2 ■ x 2 + y) = 0 and y( -► 

V 4 , - 2 = 2 / 

ans<i>*9v A 2l 


MlilN FiliP liUTD 

FDL 7/SO 


Bgjgsga i 

■ deSol','i?[2 ■ x • y - y 1 ■[2-x j;:: + yJ = 0 and y( -► 



ans<l>*9v A 2 


MlilN FiliP liUTD 

FDL 1/ 7 






■ free_off — ^ — ,x,y| 

l2-x 2 + y J 

ntfct.|-^-,y| 


2-vr + y 


intf ct <~3/v. v>l 






_ . _ , l' , 2-x-y 

■ deSo 1 ve y = = , x , y 

L 2-x 2 + y 



■■■solve<v l =2x*v/<2*x A 2 + v> ,x,v>l 

MAIN RAD AUTD FUNC S/SO 


^Ifli^SWcaYcIotherTppqklolciean UpI 

■ free_of 

x 5 + a-x 4 + x 2 y + y 3 


*-(* 2 

+ u 2 ) "”1 

-4 (x 2 + y 

2 ) 

x 5 + a . x 4 + x 2. 

-4 

X 

y + y 3 



■IAIN RAD AUTD FUN 'I ? ■ u 






1 deSolyety 1 = ode , x , y) 

-(s-x 5 + S'a'X 4 - 6-x 2 'y - 2'y 3 ) 


1 Ir 1 yebraicaTcIotheriprai', I ole 1 ea^T IJpT 1 

a x-siixyj + y cost.y. 
y ■ sin(y) - x-cos(y) 

a ode 

-(x ■ sin(y) + y ■ cos(y)) 

■ free_of(ode , x , y) 

x ■ cos(y) - y ■ sin(y) 
-(x ■ cos(y) - y • sin(y)) 

x ■ sin(y) + y ■ cos(y) 

,1 

intfctd .x>l 

MAIN RAD AUTD FUNC EVSO 


UpT ] 


■intfct(l,x) e x 

■ deSoluej y 1 = ode and y( 1) = , x , yj 

■? • y • cos(y) = <? • y ■ cos(y) - (x ■ sin(y) + y - co:^ 

■ 1 ef t( e • y • cost y) = e ■ y ■ cost y) - ( x ■ s i n( y ) 

(x • sintyj + y • costy) - sinty)) • e* = Q 

left Cans <1 >>-riffht<ans<l>>=Q 

MAIN RAD AUTD FUNC II. 


Upl 11 

2x J y z -M 

-4 x (x 2 - y 2 ) 1 
2-x 3 y - 1 
4 -(x 2 - y 2 ) ■ y 

■ t t ee_Ot 1 -r , X , '-1 

L x - 2 x 2 y^ 

. 2 x y 3 - 1 

... < 2 x A 3 y A 2 -y > / < x-2 x A 

2y A 3>,x,y>l 

IMAIN RAD AUTD FUNC ii/iO 1 


fulfil •SEr.TcSTclofiSrTfrSnlofel JST Upl 

f 2-x 3 y 2 - y 

■ deSolve y 1 = -=^ — ^ 

{ x - 2-x^ y^ 1 

x 3 - 

4 (x 2 - y 2 ) ■ y 

2 - x y 3 - 1 
x, yj 

y + x y3+ 1 = e6 

... A 3 »y A 2 -y > / < x-2 »x A 2 »y A 3 > . x . y > 

MAIN RAD AUTD FUNC iH/JO 


free_of and intfct are functions 
according to the functions 
which could be found in earlier 
Derive s ODE1.MTH. 

(download!) 


Ex 1 5 : (x - 2 y + 5) dx + (lx - y + 4) dy = 0 

a) Find the general solution using (x 0 [y 0 ) and using a parameter c. 

b) Find the special solution with v( 1) = 1. 

c) Sketch the integral curve and the direction field. 




1 deSo 1 ve y 1 = 


"(x - 2 ■ y + 5) 

2 x - y + 4 ’ 


_ ~(x - 2 ■ y + 5) 

l_ 2 ■ x - y + 4 

,e(v' = ~<x-2y+5V<2x-y+4> .x.y>l 




jRLDTS i 

tQ=l. 
xyr-^f 


■2_yl+5) 


2 ■ t - yl ■ 

y i 1 = 1 
y2'=l 
yi2= 
y3'= 
yi3= 

J&l 

w2 1 <t >= 

SELECT ONE 1ST-DRDER FUNCTIDN □NLV 


[^Iz^T^LdRe^aPhlMaiyn^VIlC 1 

^ X X X X 
V -. V -. V -. V -. V -. 

. ) 

VA \ 1 / 
s\\x- 

>-* y \ \ % 
1 \ \\\ 
1 U\N 

S S X X X X X X X --- 
v\\ X X X X X X 

x x xNt x xxxxx 

\ X X X X X X 

<***/ 1 1 

\ \ \ w 

XXxxx ■vXx x x 

X///I 

\ \ \ \ \ 

X X X X x x x 

/ / i 1 \ 

\ \ \ \ X 

X. X X XXxxx \ 

MAIN 

RAD AUTD 

DE 


i g=pr 


tnax=2G . 
tstep= . G5 
tplot=l . 
xnin= - 4. 
xnax=4. 
xscl=l . 
ynin="l . 5 
ynax=3. 
yscl=l . 
ncur i , i es=6. 
diftol-.001 
f i.-i re= . =? ri 


MAIN 


DE 


deSolve does not work. I investigated how Derive does and 
first solved the equation step by step to pack later the solving 
procedure it into one function linfrac. 

(download). 
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-(x-2-y + g) -(x-2-v + S) 

2 ■ x - y + 4 ' J - h 2 ■ x - y + 4 

1 solve(getNunCder) = 0 and get-Denom(der) ► 
x = "1 and y = 2 
den | x = x_ + ■ 1 and y = y_ + 2 

-(x_ - 2 • y_) 


2 ■ x_ 


de li- 


upT l 


f , -(x_-2y_) ) 

{ y -' = 2x_-y_ 

— ^ ^ t = ln(x_) + r 

■I k X ~ \ = ln(X-) + 


ans<l>*2 


MAIN F;HD AUTD PE fi/30 


upT 1 


ln[ " tx ~_ -3 ln[ X ~ !i + y ~ ] = 2 (ln(xJ ► 

. f in ( ~ (x ;r y-) ] - 3 ■ H 2 ^-) = 2 < 

-x_ 2 (x_-y_) _ 2.^2-®2 


(*- + y-T 


- = x_ e 


N'jminJ: Operation night introduce false solutions 


(rmi Fl* T Fl* Y FH T Y FE Y Ffi^ 

▼ 4 — - fl 1 gebra Ca 1 c Ot her Prgn 1 0 C 1 ean 

Up! | 



X- + y_) 3 



- X_ 2 " ( x_ — 

>J=X 

2 .^2 @2 



. (x_ + y-) 






x_2 






■o 

- - y-) 

.262 



(x_ 

+ y-) 3 


an s < 1 > 1 x_ 

_=x-<- 

1> and 

y_=y- 

2| 

MAIN 

FiiiP AUTD 

DE 

H/SO 



Pff hi^fcSTcloWrStioTci upI 1 

(x_ + y_) 3 



"(x - y + 3) 9 (59 


“T ■" 

(x + y- l) 3 

. "(x - y + 3) 

= f 2 '® 2 | x = 1 and y = 1 

(x + y- l) 3 

_3 = f .2.E2 

!?TTfS« M rS KTT. 1 ' : r^spj 

MAIN 

FiiiP AUTD PE i0,-'30 


fFi-Tmni Fl* Y Fl* T FH* Y FE T Fu^ 

▼ 4 — ■ fl 1 gebra Ca 1 c Other Prgn 1 0 C 1 ean 

Up! 1 


(x + y- ir 



-(x + l) 2 - 

> - y + 35 

~ i'x + :) 2 


. (x + y 

- ir 




(x+l) 2 





"(x-y + 3) 

. 2 • G8 



(x + y- l) 3 


ans <1 >/Cx+i > A 2 

MAIN 

FiAP AUTD 

PE WHf 




PffIm&aETJo&lpn? n ioI G l^ UpI ] 

(x- 

r = + - -- 1 nr 

y-ir 

■ linfrac 

(X+l) 2 

-(X 

(x 

- x + 2-y - 5 
2 ■ x - 4 ■ y + 4 ’ x ’ y 

- y + 3) 9 . ( 5 ft 

y-l ) 3 

Done 

... c ( ( " x+ 2 »v-5 > / ( 2 »x-4 v+ 4 > . x . v >1 

MAIN 

FiAP AUTD PE fi/50 



f -x + 2 y - 5 1 

LIN_FRAC , -1, 2, -5, 2, -4, 4, x, y, 1, 1 

l 2-x - 4-y + 4 J 

? 

Special case if "numerator line" and "de- 
nominator line" don't intersect. 

f -x + 2 y - 5 ) 

FUN_LIN_CCF_GEN , 2, -4, 4 

l 2-x - 4-y + 4 J 

Then substituting u = numerator leadds to a 
DE which can be treated by separation of 

x - 2 (K - 1) 3 LN(2 x - 4 y + 7) 

variables. You can follow this below. 

2 4 

f x - 2 (y - 1) 3 LN<2 x - 4-y + 7) \ 

= x + c 1 8 

l 2 4 ) 

You might extend the Tl-program linfrac 
in such a way that this special case is in- 

4 - (x - 2 (y - 1)) - 6 LN<2 x - 4-y + 7) = 8 (x + c) 

eluded instead of answering "non appli- 

- 6 LN(2 x -4y +7> +4x-8y ♦ 8 =8x +8c 

cable". 

(- 6 LN(2 x - 4 y + 7) + 4 x - 8 y + 8 = 8 x + 8c) - 8 x 


The final solution: 

(download) 

- 6 LN(2 x - 4-y + 7) - 4 x - 8 y + 8 = 8c 

Derive provides a special utility function 
for this case. 


Ifi 1 ftebcalcaTtlotherlprqm I ole: 1 ean UpI 1 

■ solve(u 
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■ deSolve 

^ u + x + 5 

- _ x + 2-y 5,y) y “ 2 

-x + 2-y-5 . u + x + 5 

2 ■ x - 4 ■ y + 4 l y 2 

u' + 1 -u 

2 

u' + 1 "u .. ' 

2 ■ ( u + 3) 

2 - 2-( U + 3)’-'’ U J 


MAIN 

RAP AUTD PE B/10 
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3 1 ln(2-u + 5) _u _ _ 

4 + 2 - Li - x 

, 3- ln(2 u + 3) . u i 

t +- 2 =liil-x|u = -x + 2 ■ y - 5 


3- ln( -(2-x - 4-y + 7)) 


- 2-y + 5 . 


gEgggjg up] ] 


3- ln( -(2-x - 4-y + 7.1) x - 2 ■ y + 5~ 


. < r ,i;-(2-x-4-y + 7» _ x-2 2 y + 5 =el _ x j. s 

2-(3 - ln( "(2-x - 4-y + 7)) -2-(x-2-y + 5)) > 
■ ejtpand(2 (3 ■ ln( -(2 ■ x - 4 ■ y + ?)) - 2 (x - 2 ► 
4 2 - x - 4 y + 7)) - 4 - x + 9 y - 2Q = 3 6 1 - 3 x 

■liilH FiiiP AUTD PE B.'jQ 


p 
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iSKErass upi ii 

— ■ ■ J 1 1 II, l, i 1 a -r 1 y I i j) J. L.A 1 — i y I -Jf) - 


( J 1 1 n, { i- a t‘ y ' i f) J- i, a‘ — j. y ' -Jj) 

expand(2 (3 ln( -(2 ■ x - 4 y + 7)) - 2 (x - 2 ► 


42x-4y + 7))-4x + 8y-20 = 8ei-8x 
■ (6 1 n( ■( 2 ■ x - 4 - y + 7)) - 4-x + 8-y-20 = 8| 
6- InC -(2-x- 4-y + 7® + 4-x + 8-y - 20 = S i> 


■(0- ln( -(2 x - 4 y + 7)) + 4 x + 8-y - 20 = 8> 
6 ■ 1 nf -( 2 ■ x - 4 ■ y + 7)) + 4 x + S y = S Sl + :► 

...<2*x-4*w+7>>+4*x+8*w=8*Pl+2Q 

MHIH FiiiD HUTD PE 


: 1 infr-ac{r,u,'wO 


: Prgpi 

s Local n,d,s,c,f ,cc,J0,r_, 1_ 

: get-Nun < r "? -*n : getDenon < r !> -*d 

o?-d<n,o)*d<Id,uI>-*s 

s If s=G Then 

:Disp "non applicable" 

:0oto end 
: End I f 

: ~nlu=0 and o=0+c 
: ~dlu=0 and o=0*f 
: <d<n,u)*f-c*d<d,u))/s+JE1 


Ex 16: Clairaut equation y = xy' + y’ - y' 2 

Find the general solution, give the singular solution and plot the respective graph. 
Find the special solution(s) through R(2|l). 


^IhlYu,, ',nn77, -7: i, UP i n 

■ clair(x ■ o - 

y = y 2 -y) 

"GenSol" cx-y-c 2 + c = 0 
"flux" -2o + x+ l= G 

(x + l) 2 

"SingSol" y = ' 4 J 

c 1 a.i r < 

-y=u A 2-u> 

1 1 

IMHIH FiHD HUTD FUH'I l.^>> [ 


f»gl R ,.aLx 3 ic^:!o^rlp f -a,IQ]ci,gT UpT i 


"c^ + 3-c - 1 


1 zeros( _ c 2 + 3 ■ c - 1 , c) 


-(J5 ~ 3) J5 


H 


2 2 

■ x-x-x2 + r .| r . = U^jl Jg±3^ u ^ 

< J<5 >-3 V2 . < J<5 ?*3 >/2>*v3 <x?l 




J-LDTS 

Plot 1 :Lil nr-cir-a 

. (x + l) 2 
*V 1 = 4 

xy2=c -x-c 2 + c|c = sey<; k , k , ■ 

2 l ( "(is - 3) 
yy3=xc - c A + c | c = | ^ 

y4=* 

v4<x>= 


^ 1 

,3,.5) 

J5 + 5 j 



: Local p,q,pi ,p2, sing, sings 
: lefttr J-»-p:rignt<r>-»-q 


di-nnjif FT? TfrS^ 

| ▼ E— ■ JContr o 1 1 1 /Q |Uar |F i nd... |Mod< 
s ■: -lairtrOI 
:Func 

: Local p,q,pl 
t<iO-*p:ri 
=0Io=c-*p1 

:x*l irnit<d<l iDiit<p,y,x*o-y),y),y,x*o-y)- 
dtq,o)=0*p2 

:zeros< left(p2> ,v) l 1 1+sing 
: limit(r,o,sing}+si ngs 

: [ [ M 0enSol r, ,pl] [ "flux ff ,p21 [ "SingSol" ,solo 
e ( s i ngs , y ) 1 1 


deSolve cannot be applied. 
Flelp yourself and write a 
function clair. 

(Replace y' by v.) 


Ex 17: Which is the curve with its tangent's segment between the axes having constant length a = 21 
Give both general and singular solutions. 


f^EralCal 


^ Y FH- Y FE Y Fb T Y \ 
1 c |0ther- 1 Ft- g n IOjCl ean Up | | 


"GenSol " 
"flux" 


27T „ 

c-x - y - i „ = 0 


Jc 2 H 


, 2 + l) 


X 2X3 _ 2 2^3 

1-llr.Qtr.l < n amrl 


clair <x*a-v=2aAT <l+a A 2 > > 


i^IfllSnaKYclo&lrn^IoIcie^ UpT 1 

L 




► 

3 _ 2 2/3 ! 


2 2X3 _ X 2X3 


f 

, 2 x 3 1 0 ■="" J x 


x 1 '' 3 


1 


uhen 

r x 2/3 _ 2 2/3 

< 0 


x 2 -- 3 

ans <1 > [ 3 ,2 ]| 

MHIH FiHD HUTD 


FUH'I i.’OO 




f^glHliibraT^R^ 


PrgrclOlCleai 


iTu^in 


L 2 + l) 


U 2 -!) 


- = 0 | o = tan(t) 


X - 2 ■ |(oos«iJl 3 | = 0 


V^IhiY-U,, 'Y in .!,:.. - 7 : ii:l..'i';T Hi- i 


- 1 y = tanCi^ and x = 2-fico^. 


2 ■ sin(t) ■ cosCt) • |cosCt)| - y = 2 ■ sin(t) ■ si<^ 


■ soloe(2 ■ sinCt) cos(t) ■ |cos(t)| - y = 2-sin^ 
y = 2 ■ sin(t) ■ cos(t) ■ |cos(t)| - 2 ■ sin(t) ■ si<^ 


■ DrawParn 2 • If cost! till 2 s i nC t) • cos( t) • 1 

... <t >*si?rn<cQs<t >>. ~n/2.n/2 . .1 

MHIH RHP HUTD FUHC fi.'jp 


[▼ Ifi 1 yebr^afcaYc [other- Ipry'n I □ III: 1 ealT Up! 1 

■ 

Cla + U - y = Jl + « 2 ] 

" GenSo 1" c ■ x - y + -j= 

Jc 

2 

"flux" 

(y 2 + i) 3 ' 

► 

J- + X = 0 

clair<x*u-v=-2*u/<J'<i + u A 2 >>>l 

MHIH 

FiHD HUTD 

FUH'I F75':' 


PiSIffiKitIT-KIi hiZAl 

,J 1 

j-ldt; 

Plot l:LL □ x : d h : c£ 


•^y 1 =c ■ x — , = | c = seq( k , k , 

Jc 2 + 1 

- 4 , 4 , 1 ) 

■ y 'y2=c ■ x + ■ | c = seq( k , k ? 

Jc 2 + 1 

y 3 =i 

y 4 = 

y 5 = 

- 4 , 4 , 1 ) 


v3<x>= 

MHIH RHP HUTD FUHC 



More explanations can be 
found in DNL#3. 




PiSluSraTtaiofflSrWS.IoTciegS' UpI 


"SingSoll 11 

y = 1 1 ±: 1 + 1 1 ■ x + G1 2 1 


"SingSol2 M 

G1 + 1 = G1 


"GenSolx" 

x =-%- - 2 i + 2 
e - 


"GenSoly" 

y = c 'V c -t 2 + 2 
e ^ 

secr<k.k. 

"5.5. .5>-H; 

MAIN 

RAD litlTD 

FAR ££^30 


[^ImSraKTcIn&lPrlfnloIci^ UpI 1 

.IW” 

+ c 


Xp_dx_ 

f auxlCp-j q_j x_, 



Done 

■ lagrfy = x 

+ x-y + y 2 ) 



"SingSoll" 

y = (El + l)x + El 2 


"SingSol2" 

ei + 1 = ei 

... <r <b_. > 

c_>> Aauxl <p_ . _ x_ . v_ > 

MAIN 

RAD liUTD 

FUNC ££/30 




lagr(equ) with v for y' follows the Derive function from 
ODE1.MTH. You find the explanation in DNL#3. 

(download) 

No singular solution, a family of special soultions. 




DelUar o 
1 agr(y = v^-3-y + x) 

"SingSoll " y = x + @l 2 
"Sing5ol2" 1 = G1 


"GenSolx 


3 ■ 


la?rr<v=u A 3-3*u + x> 


[▼ ^51fh oetoralca lclotherlprSi I ole 1 ean 

Upf 


3 ■ t 2 

"GenSolx" x= — — + 3 

t + C 


-T 3-^ 




+ c 

■ DrawFunc 

- 2 

□one 

■ seg( k , k , _ 3 , 3 , . 5) + c 
€-3 -2.5Q -2.00 -1.5G -1.00 

-.5* 

se<r <k . k . 

-3.3. .5>+c 


MAIN 

RAD hUTU FAR £6^30 





Singular solution v = x - 2 together with a family of special 
solutions. 


f^k&aETc 

# FH^ Y FE Y FfiT Y ' 

Other iPngrn 1 0|C lean Up| 


"SingSol 1 " y= JfI _ 1 
"SingSol 2 " ( = El 

"OenSolx" x = c - -j=- 

1 asrr<v=u#x/C 

Xu>-1>> 


MAIN RAD AUTI 

FUh 

■I £B/30 



1 7 |r 1 gobralc ale lothor-lpr qn ! o(c 1 aan UpI I 

■ solye((x ■ v + y) 2 = y 2 
v (Jy + l) • 

y f y) 

-y ■ (Jy - l) ■ x 

1 w 1 

■ 1 agr y = — 

v - 1 Wl 3 v - 1 

u- 1 J 

...-. . .. E1-CJU+ l)-xl 

SingSol 1 y m1 _ , | 


■IAIN 

RAD AUTD 

FUNC £B/3D 



[1] Giinter-Kusmin, Aufgabensammlung zur Hoheren Mathematik, Berlin 1964 

[2] The Calculus Problem Solver, Research and Education Association 1985 

[3] Biiktas, Aufgabensammlung, Diesterweg 1978 

[4] Gunter Redl, T A 3 -Course Materials 
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Cellulara automater / Rule 30 

I leave David's part in Swedish, because I beleive that it is more or less self explanatoiy. Some addi- 
tional comments are in English (red). 


En cellular automat ar en samling “fargade” celler i ett rutnat. Fargema uppdateras enligt en speciell 
regel. 

Vi tanker oss en rad, med obegransad utstrackning at saval vanster som hoger, med rutor eller celler, 
som ar fyllda med antingen svart eller vit farg. 



Vi tanker oss aven att varje cell i raden andrar farg enligt en regel beroende pa dess och de narmaste 
grannamas farg. Fargandringen sker samtidigt i alia celler. Regler for fargandring kan anges pa en del 
olika satt. Pa detta satt ger regeln upphov det olika generationer. 

Vi ko mm er att visa att en undersokning av nagra cellulara automater med t.ex. DERIVE, Excel, Voyage 
200, Mathematica and Wolfram Research's A New Kind of Science Explorer utgor ett intressant pro- 
jektarbete inom ramen for kursen Matematik Diskret vad galler bl.a. binara tal och satslogik. 

Josef Bohm ar ordforande i Derive User Group och har i manga ar verkat som gymnasielarare i Oster- 
rike. 

David Sjostrand, lektor i matematik vid Elof Lindalvs gymnasium, Kungsbacka. 


En kand och intressant regel ar rule 30, som kan defmieras sa har: 


ride 30 


■H" 

■H D 

m 

■ 



■ 

■ 

■ 

s 


4 ^ 


00011110 


Ovning. Det finns 256 regler av denna typ. Visa detta. 

Pa http://mathworld.wolffam.com/ElementaryCellularAutomaton.html finns alia de 256 elementara 
endimensionella cellulara automatema askadliggjorda. 

Ovning. Bestm antalet funktioner f : A B , dar antalet element i A ar a och antalet element i B ar 


b. 

Om man identifierar vitt med 0 och svart med 1 kan man formulera en regel av denna typ med en 
funktion / : {0, 1} x {0, 1} x {0,1} — > {0, 1} pa foljande satt: 


Fargen pa en cell bestams av f(p,q,r ) dar p, q, och r fargen pa den vanstra grannen, cellen sjalv re- 


spektive den hogra grannen. For rule 30 galler saledes att 

/(l, 1 , 1) = 0 , /( 1 , 1 , 0) = 0 , /( 1 , 0 , 1) = 0 , /( 1 , 0 , 0) = 1 , /(0, 1 , 1) = 1 , /( 0 , 1 , 0) = 1 , /( 0 , 0 , 1) = 1 , /(0, 0 , 0) = 0 


Vi ser att (1, 1, 1), (1, 1, 0), (1, 0, 0), ... (0,0,1), (0, 0, 0) svarar mot talen 7, 6, 5, ... 1, 0 i binar form. 
Vi ser aven att hogerleden bildar talet 00011 110, vilket ar 30 i binar form, darav namnet rule 30. 
Man kan lyckligtvis hesta mm a ett polynom i p, q och r, som ger vardena av f. 


A. Ansatt ett polynom 

Ett satt att gora detta ar att ansatta 
/ (p, q, r) = apqr + bpq + cqr + drp + ep + fq + gr + h 

samt att sedan bestamma de 8 koefficientema , a, b, ..., h genom att losa det ekvationssystem som 
bestar av de 8 likhetema ovan. Man kan med fordel anvanda DERIVE. 
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InputMode := Character 


ff(p, q, r) := a-p-q-r + bpq + cqr + drp + ep + fq + gr + h 

S0LUE( [ff ( 1 , 1, 1) = 0, ff(1, 1, 0) = 0, ff(1, 0, 1) = 0, ff(1, 0, 0) = 1, ff(0, 1,1)= 
1, ff(0, 1, 0) = 1, ff(0, 0, 1) = 1, ff(0, 0, 0) = 0], [a, b, c, d, e, f, g, h]) 


[a=2^b=-2^c=-1^d=-2^e=lAf=1 A g=1^h=0] 

ff30(p, q, r) := 2*p*q*r - 2xpxq - q*r - 2*rxp + p + q + r 

Vi ser att / ( p , q,r)-2 pqr -2 pq - qr - 2 rp + p + q + r . 


We can try another approach, interpreting the triple p,q,r as a binary number pqr, taking its decimal 
equivalence and finding the respective polynomial dependent on only one variable: 

AV = 0,/6) = 0,/5) = 0/4) = 1/3) = 1/2) = 1/1) = 1/0) = 0 


7 G 5 4 3 2 

fb(x) := a x + b'X + c x + d x + e x + f-x + g x + h 

S0LUE( [ fb( 7) = 0, fb(G) = 0, fb(5) = 0, fb(4) = 1, fb(3) = 1, fb(2) = 1, fb(1) = 1, fb(0) 
= 0], [a, b, c, d, e, f, g, h]) 


1 23 41 89 481 1787 7 

a - - s, b — c “ “ ^ d — /\ e “ “ ^ t ^ g — ~ 

3G0 3G0 72 3G 90 360 12 


7 

x 23 G 41 5 

f b30( x) := - + xx - xx + 

3G0 3G0 72 


89 4 481 3 1787 2 7 

xx - XX + XX - XX 

36 90 360 12 


7 0 
6 0 
5 0 


TABLE(fb30(x) , x, 7, 0, -1 ) 


4 1 
3 1 


2 1 


1 1 


0 


0 
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B. Logik 

Om vi identifierar svart med SANT och vitt med FALSKT ser vi att vi kan uttrycka rule 30 pa detta 
satt 

[p a -i q a — i r ) v (— i p a q a r) v (— 1/> a q a — ir) v [—^p a — i<? a r ) . (*) 

Om vi gor foljande definitioner i DERIVE, kan vi bestamma ett polynom, som ger rule 30. 

and_(p, q) := pq 

°r_(p, q) := p + q - and_(p, q) 

not_(p) := 1 - p 

or (p , q, r, s) := or_(p, or_(q, or_(r, s))) 

and (p, q, r) := and_(p, and_(q, r)) 

Now we translate Rule 30 from above (*): 

or (and (p, not_(q), not_(r)), and (not_(p), q, r) , and (not_(p), q, not_(r)), 

and (not_(p), not_(q), r)) 

4 2 2 2 2 3 2 

p q r (q - 1) (r - 1) + p q r(1 - r)(q - 1) (3 q r(r - 1) + qr(5 - 3r) - 1) + 

242 232 2432 

p ' (3 • q r (r - 1) + q r (1 - r)(6r - 11) + q -(3r - 1 1 ■ r + 6r + 4r - 1) + 

2 42232 24 

q ■ (4 ■ r - 6 ■ r + 1) - r(r - 1)) - p • ( q r (r - 1) + 2q r (1 - r)(r - 3) + q (r - 

32 2 2 32 22 

6 ■ r + 5 ■ r + 3 r - 1) + q (3 r - 8 r + 3) - r + 3 r - 1 ) + q ■ r • (r - 1 ) + q r (2 - 

r) + q (1 - 2 r) + r 

After expanding we find a very bulky expression, which can be "simplified", because alle powers of p, 
q and r are 1. We use the asterisk as TimesOperator, because we intend to transfer this formula to 
Excel. 

444 443 442 434 433 432 424 

p q r - 2p q r + p q r - 2p q r + 4p q r - 2p q r + p q r - 

423 422 344 343 342 334 333 

2- p q r + p q r - 3 p q r + 6 p q r - 3 p q r + G p q r - 14 p q r 

332 324 323 322 32 3 2 3 

+ 8 ■ p q r - 3p q r + 8p q r - 4p q r - p q r - p qr + p qr + 

244 243 242 234 233 232 

3- p q r - G p q r + 3 p q r - G p q r + 17 p q r - li p q r + 

224 223 222 22 22 22 2 2 

3p q r - lip q r + 6p q r + 4p q r - p q + 4p qr - Gp qr + p q 

222 44 43 42 34 33 32 

- p r + p r - pq r + 2pq r - pq r + 2pq r - 8pq r + Gpq r - 

24232222 2 2 

pq r + Gpq r - 5pq r - 3pq r + pq - 3pqr + 8pqr - 3pq + pr - 

33 32 23 22 

3pr + p + q r - q r - q r + 2q r - 2qr + q + r 
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pqr - 2pqr + pqr - 2pqr + 4pqr - 2pqr + pqr - 2pqr + pqr - 3pqr 
+ 6pqr - 3pqr + Gpq-r - 14pqr + 8pqr - 3pqr + 8pqr - 4pqr - pqr 

- pqr + pqr + 3pqr - Gpqr + 3pqr - Gpqr + 17pqr - 11-pqr + 3pqr 

- 11-pqr + 6-p-qr + 4pqr - pq + 4pqr - Gpqr + pq - pr + pr - pqr + 
2pqr - pqr + 2pqr - 8pqr + Gpqr - pqr + Gpqr - 5pqr - 3pqr + 

p q - 3-p-q-r + 8 p q r - 3pq + p r - 3 p r + p + qr - qr - qr + 2 q r - 2 q r + q 
+ r 


TimesOperator := Asterisk 


p*(q*(2*r - 2) - 2*r + 1) + q*(1 - r) + r 

Compare with expression ff30 from above 

f f 30( p , q, r) = p^(qx(2xr - 2) - 2 *r + 1) + qx(l - r) + r 
rul©30(p, q, r) := 2*pxq*r _ 2xpxq - 2xpxr + p - q*r + q + r 

Denna formel ar onodigt stor. Man gora ersatta den med en mycket enklare genom att utnyttja att 

n n n 

p = p, q = q,r =r,n>\. 

Ovning. Enligt S. Wolfram, A New Kind of Science sid 884 ges rule 30 av 

^ ^ ) . v star for XOR (exlusive or). 

a) Visa att 

( p A—<q a — ir) v (-1 p a q a r) v (— 1/7 a q a — 1 r) v ( — 1/7 a — 1^ a /*) <=> pv(q v r) genom 
att anvanda funktionen TRUTH TABLE i DERIVE. 

b) Anvand uttrycket /;v (<;/ v r ) for att skapa ett polynom for rule 30. Anvand DERIVE eller 
rakna for hand. 


exprl := (p a 

J 

_Q 

> 

J 

1 

< 

(- , P*q*r) 

- (-■ P 

> 

_Q 

> 

J 

r) v (-1 p ^ 

-> q ^ r) 

expr2 := p XOR 

(q ^ r) 









P 

q 

r exprl 

IFF expr2 




true 

true 

true 

true 




true 

true 

false 

true 




true 

false 

true 

true 

TRUTH_TABLE(p, 

q, r, exprl 

IFF expr2) = 

true 

false 

false 

true 




false 

true 

true 

true 




false 

true 

false 

true 




false 

false 

true 

true 


false 


false 


false 


true 


p 
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C. Rule 30 med Excel 

Vi kan nu askadliggora ett antal generationer av Rule 30 i ett kal- 
kylblad i Excel. I den fiorsta generationen Finns endast en svart cell 
eller en cell innehallande vardet 1. 





Ml 






1 

Ml 

Ml 




i 

1 

0 

0 

hr 



i i 

0 

Ml 

M 

i i 


hi 

1 0 

0 

Ml 

0 

0 0 



Anvisningar 

1 . Skapa en egendefmierad funktion rule30Q?, q,r) = 2 pqr -2pq - qr - 2 rp + p + q + r . 

a) Oppna Visual Basic Editor genom att valja Verktyg Makro Visual Basic Editor. 

b) Valj sedan Inf oga Modul. 

c) Skriv in foljande 

Function rule30 (p, q, r) 

rule30 =p+q+r-q*r-2*p*q-2*p*r+2*p*q*r 
End Function 

Funktionen Rule30 kan nu anvandas i ett kalkylblad. 

2. Skriv in 1 i cell DG1 samt =Rule30 (DF1 ; DG1 ; DH1 ) i cell DG2. Kopiera sedan formler sa att du 
far ett kalkylblad av den typ som fins i figuren ovan. 

3. Du kan sedan modifiera cellemas storlek med Format Kolumn Bredd, Format Rad Hojd. 
Man kan latt man gora de celler, som innehaller 1 :or svarta. 

I figuren har bredvid har vi askadliggjort 78 generationer av Rule 30. 



You can format the cells de- 
pendend on their contents. 

All cells with content 1 are red 
(font and pattern), all cells with 
content 0 are full white. 
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Ovning. Askadligor ett antal generationer av rule 90 och rule 110 med Excel. 

ff(p, q, r) := a*p*q*r + bxp*q + cxqxr + dxrxp + exp + f*q + gxr + h 

S0LUE( [ f f ( 1 , 1, 1) = 0, ff(1, 1 , 0) = 1 , ff(1, 0, 1) = 0, ff ( 1 . 0, 0) = 1, ff(0, 1,1)= 
1, ff(0, 1, 0) = 0, ff(0, 0, 1) = 1, ff(0, 0, 0) = 0], [a, b, c, d, e, f, g, h]) 


[a=0^b=0^c=0~d=-2~e 


f = 0 ^ g = 1 


^ h = 0] 


- 2xrxp + p + r 



This are Excel-creations 



Med hjalp av DERIVE kan man exempelvis astadkomma dessa figurer. Den ursprungliga generationen 
bestar av flera l:or och 0:or. 


A New Kind of Science Explorer 

A New Kind of Science Explorer ar en programvara fran Wolfram Research med vilken man kan gora 
alia de experiment med cellulara automater vilka beskrivs i [ 1 .] Man kan “zooma” var som heist i en 
cellular automat. 



Med programmet kan man aven undersoka cellulara automater av hogre dimension och med mer kom- 
plicerade uppdateringsregler an vad vi beskrivit ovan samt sadana som ar uppbyggda av fler farger. 

Litteratur: 

[1.] S. Wolfram, A New Kind of Science 

Lankar: 

http://mathworld.wolfram.com/ElementaryCellularAutomaton.html 
http://www.math.usf.edu/~eclark/ANKOS reviews.html 
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D. A possible Derive implementation: 


binar(n, n_ , i_, j_) := 

Prog 

OutputBase := 2 
n_ := STRING(n) 

OutputBase := 10 

i_ := 1 

#1: j_ := 8 - DIH(n_) 

Loop 

If i_ > j_ 

RETURN VECTOR (k_ - 48, k_, NAME_TO_CODES(n_)) 
u_ := APPEIID("0" , n_) 
i_ :+ 1 


#2: [binar(30) p binar(llO), binar(90)] 


00011110 

01101110 

01011010 


cell_aut(r, u , st_, stt_, i_, k_, 1_, pts_, r_, ant) := 
Prog 
i_ := 1 

pts_ := [[0, 100]] 

1 _ := -1 
stt_ := [1] 

Loop 

st_ := APPEHD([0, 0], stt_, [0, 0]) 
stt_ := [] 
k_ := 1 

#3 : Loop 

r_ := st_ik_-4 + st_|(k_ + l)-2 + st_i(k_ + 2) 

If (binar(r)) r(8 - r_) = 1 

pts := APPEND(pts_, [ [ k_ + 1_ - 1, 100 - i_]]) 

stt := APPEIID(stt_, [ (binar(r)) i(8 - r_)]) 

k_ :+ 1 

If k_ > DIM(st_) - 2 exit 
If i_ = n 

RETURN pts 

1_ 1 
i_ :+ 1 


The initial number is 1. For finding the next gen- 
eration I append a double zero in front and at the 
end of each sequence (in red). 

Plotting starts at (0| 1 00) and moves downwards - 
like a curtain. 

All following generalisations base on 

cell aut (rule#, #ofGenerations) . 


0 0 10 0 
0 0 1110 0 
0 0 1 1 0 0 1 0 0 
0 0 110 11110 0 
0011001000100 
001101111011100 
011001000010010 
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#4: cell_aut(30, 120) 

#5: cell_aut(90, 120) 

#6: cell_aut(110, 120) 




#7: cel l_aut(3i, 60) 
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In my first generalisation I start with an initial binary sequence - self written or produced randomly. 


VECTOR (RANDOH(2), k, 100) 

[ 1 , 1 , 1 , 0 , 0 , 1 , 1 , 1 , 1 , 0 , 1 , 1 , 1 , 0 , 1 , 1 , 1 , 1 , 1 , 0 , 1 , 0 , 1 , 1 , 0 , 0 , 1 , 1 , 0 , 1 , 


o, 

1, 

0, 

1, 

1, 

1, 

0, 

1, 

1, 

o, 

1, 

1, 

o, 

i, 

o, 

o, 

1, 

i, 

o, 

o, 

o, 

0, 

o, 

1, 

i, 

0, 

1, 

1, 

1, 

0, 

1, 

0, 

1. 

1, 

o, 

1, 

o, 

1, 

1, 

1, 

0, 

0, 

1, 

1. 

1, 

1, 

1. 

1. 

1, 

o, 

0, 

1. 

1, 

0, 

0, 

0, 

o, 

1, 

1, 

1, 

o, 

o, 

1. 

o, 

o, 

1, 

o, 

1, 

1, 

1] 


















testl := 

[1. 

1, 

1, 

0, 

0, 

1, 

1, 

1, 

1, 

o, 

1, 

1, 

1, 

0, 

1, 

1, 

1, 

1, 

1. 

o, 

1, 

0, 

1, 

1, 

0 

, 0, 

1 

1, 

o, 

1, 

0, 

1. 

o, 

1, 

1, 

1, 

o, 

1. 

1, 

o, 

1. 

1, 

0, 

1, 

o, 

0, 

1, 

1, 

0, 

0, 

o, 

o, 

0, 

1, 

1, 

0, 

1, 

1, 

1, 

o, 

1, 

o, 

1, 

1, 

o, 

1, 

o, 

1, 

1, 

1, 

o, 

o, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

o, 

0, 

1, 

1, 

o, 

o, 

0, 

o, 

i, 

1. 

1, 

o, 

0, 

1, 

o, 

o, 

1, 

o, 

1, 

i, 

1] 
















cen_aut_gen(30, testl, 100) 


cell_aut_gen(37 , testl, 100) 



test3 and test4 consist both of 99 zeros and one 1 on different places. See the difference applying 
Rule90 (next two pictures). 
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test3 := 

[0, 

o, 

0, 

o, 

o, 

o, 

o, 

0, 

o, 

o, 

0, 

o, 

o, 

0, 

o, 

o, 

0, 

o, 

o, 

0, 

o, 

o, 

0. 

o, 

o, 

0. 

o, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

o, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

o, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

o, 

o, 

0, 

o, 

o, 

0, 

o, 

o, 

0, 

o, 

o, 

0, 

o, 

o, 

0, 

o, 

0, 

0, 

o, 

0, 

0, 

o, 

0, 

0, 

o, 

0, 

0, 

0, 

0, 

o, 

o, 

0, 

o, 

o, 

0, 

o, 

o, 

1, 

o, 

o, 

o, 

0] 















test! := 

[0, 

o, 

0, 

o, 

o, 

0, 

o, 

o, 

0. 

o, 

0, 

o, 

o, 

0, 

o, 

o, 

0, 

o, 

o, 

0, 

o, 

o, 

0, 

0, 

o, 

0, 

o, 

0, 

0, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

o, 

0, 

o, 

0, 

0, 

0, 

o, 

0, 

0. 

o, 

0, 

0, 

o, 

0, 

0, 

o, 

0, 

0, 

o, 

0, 

o, 

0, 

0, 

o, 

0, 

0, 

1, 

0, 

0, 

o, 

0, 

0, 

0, 

0, 

o, 

o, 

0, 

o, 

o, 

0, 

o, 

o, 

0, 

o, 

o, 

0, 

0] 

















cell_aut_gen(50, test3, 100) 



cell_aut_gen(122 , test3 , 100) 



In my last generalisation of this kind of cellular automates I include generating the first generation 
randomly in the program. cell_aut_rnd (rule#, #ofColumns, #ofGenerations) produces sur- 
prising results. 
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Rules 90, 80 and 137 



Try special rules and observe the various kinds of the created patterns: 

rules 8, 136 
rules 4, 37, 56, 73 
rules 18, 45, 146 
rule 110 


Go on and design your own wallpaper. 
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One example to convert from DERIVE6 to DERIVE5 

(No responsibility for possible mistakes can be taken.) 


The Derive 6 mth-file trrules . mth 
(containing some greek letters!) 



#4: Cosine Pule, given a,b and y_ 

2 2 2 

#5: c = a + b - 2-a-b-C0S(y_) 

as it appears in Derive 5 



You can see that the greek letters are replaced by strange other characters - cuased by using Unicode 
in Derive 6. 


I have d6tod5 . exe in the same folder as trrules .mth and run d6tod5 trrules: 
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F:\D6tod5.exe 


File TRRULES5 . MTH created 
Press any key_ 


I could also open the DOS-Window and then call the program. 


The next key press should terminate the little step into DOS past. 


Eingabeaufforderung 

Microsoft Windows XP [Uersion 5.1.2600] 
(C) Copyright 1985-2001 Microsoft Corp. 

C:\Dokumente und Einstellungen\User>f : 

F:\>d6tod5 trrules 

File TRRULES5 . MTH created 
Press any key 

F : \>exit 


I leave DOS by typing exit. 

We can see that TRRULES 5 .MTH has been created. Opening this file with Derive 5 shows the con- 
verted file, containing the correct greek letters. 
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Albert Rich 

Hello Valeriu, 

Thank you for pointing out the problem with Derive's trig functions when approximating large magni- 
tude arguments. I have resolved the problem as you suggested by temporarily increasing the preci- 
sion while reducing the argument modulo 2*pi. Now SIN(1 0 A 1 1 ) approximates to 0.9286936604 at 10 
digits of precision. 

Since 10 A 1 1+0.5 justifiably approximates to 10 A 11 at 10 digits of precision, both SIN(10 A 1 1+0.5) and 
your sin (1 0 A 1 1+0.5) approximate to 0.9286936604. However, in Mixed precision mode, 
SIN(10 A 1 1+0.5) correctly approximates to 0.9927992642. In fact, this is the best application of Mixed 
precision mode that I have ever encountered! 

I will answer the other issues raised in your email in a subsequent email. 

Aloha, 

Albert 


Hello Valeriu, 

I have implemented the 3 Euler Substitutions used to rationalize integrands of the form 
F(x,SQRT(a+b*x+c*x A 2)) as given by formulas in Section 2.251 of Gradshteyn&Ryzhik's "Table of 
Integrals, Series, and Products", Fifth Edition. 

Derive is now able to integrate 1/(x*SQRT(x A 2+x+1)+2*x+2) to the sum of two logs and an arctan- 
gent. Mathematica's on-line integrator gives an inexplicable antiderivative of this expression in terms 
of Rootsum. Maple V, version 4 is unable to integrate it. I would appreciate any readers of this email 
having a more recent version of Maple, try integrating this expression using Maple. 

Used indescriminantly the Euler Substitutions results in high degree integrands that give more compli- 
cated results or even hang Derive. This is because partial fraction expansion of the original integrand 
is easier than expansion of the rationalized integrand. Therefore, as a hack, I have restricted the use 
of the Euler substitutions to integrands that are pure reciprocals and b is non-zero. In an effort to nar- 
row the scope of this restriction, I would appreciate any examples that: 

1 . The Euler substitutions can get; 

2. Derive 6.01 cannot get; and 

3. are pure reciprocals or have a b=0. 

Finally, is there a better substitution to rationalize square roots of binomials of the form 
F(x,SQRT(a+b*x A 2)) or even better F(x,SQRT(a+b*x A n)) where n is a positive integer? 

Thanks again for your help in furthering the education of Derive. 

Aloha, 

Albert 
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Dear Derive and Tl-Community 


It is a really difficult problem. 

It does not have solution? 

We have to put the numbers 1,2,3, 19 in 

the diagram. 



The following sums have to be 38. In the verticals, diagonals shown. 





As you see, we need 

Three sums of five summands equal to 38 
Six sums of three summands equal to 38 
Four sums of four summand equal to 38. 


Milton Lesmes, Bogota, Colombia, miltonlesmes@hotmail.com 


